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^ : ABSTRACT 

^ I I show that the relativistic winds of newly born magnetars (neutron stars with 

I petagauss surface magnetic fields) with initial spin rates close to the centrifugal 

^ I breakup limit, occurring in all normal galaxies with massive star formation, can 

OO I provide a source of ultrarelativistic light ions with an injection spectrum, 

I steepening to at higher energies, with an upper cutoff at 10^^ — 10^^ eV. In- 

^ I teractions with the CMB yield a spectrum at the Earth which compares favorably 

with the spectrum of Ultra-High Energy Cosmic Rays (UHECR) observed at en- 
^ _ ergies up to a few x 10^° eV. The fit to the observations suggests that ~ 5 — 10% of 

I the magnetars are born with rotation rates and voltages sufficiently high to allow 

^ I the acceleration of the UHECR. The form the spectrum incident on the Earth 

I takes depends sensitively on the mechanism and the magnitude of gravitational 

I wave losses during the early spindown of these neutron stars - pure electromag- 

^ I netic spindown (the injection spectrum) yields a GZK feature (a flattening 

■ ■ ■ ' of the E^J[E) spectrum) below lO^'' eV, rather than a cutoff, while a moderate 

GZK cutoff appears if gravitational wave losses are strong enough to steepen the 
injection spectrum above 10^° eV. The flux above 10^° eV comes from magnetars 
in relatively nearby galaxies {D < 50 Mpc.) 

I outline the probable physics of acceleration of such particles in a magnetar's 
wind - it is a form of "surf-riding" in the approximately force free fields of the 
wind. I also show how the high energy particles can escape with small energy 
losses from the magnetars' natal supernovae. In particular, I show that the elec- 
tromagnetic energy emitted by the magnetar "shreds" the supernova envelope in 
times short enough to allow most of the relativistic energy to escape largely un- 
mimpeded into the surrounding interstellar medium, where it drives a relativistic 
blast wave that expands to parsec scale before slowing down to nonrelativistic 
speeds. I also show that since the ions are accelerated in a region where the mag- 
netic field has the structure of a strong electromagnetic wave but propagate at 



also Department of Physics and Theoretical Astrophysics Center 



- 2 - 



larger radii through a region of weaker magnetic field near the rotational equator 
of the outflow, the ultrahigh energy particles escape with negligible adiabatic and 
radiation losses. 

The requirement that the magnetars' relativistic winds not overproduce inter- 
stellar supershells and unusually large supernova remnants suggests that most of 
the initial spindown energy is radiated in khz gravitational waves for several hours 
after each supernova. For typical distances to events which contribute to £■ > 100 
EeV air showers, the model predicts gravitational wave strains ~ 3 x 10~^^. Such 
bursts of gravitational radiation should correlate with bursts of ultra-high energy 
particles. The Auger experiment should see bursts of particles with energy above 
100 EeV every few years. 

Subject headings: acceleration of particles — cosmic rays — relativity — stars: 
neutron — supernovae: general — supernova remnants 

1. Introduction 

I study the possibility that relativistic winds from rapidly rotating magnetars, neutron 
stars with surface dipolc fields on the order of 10^^ Gauss (Duncan and Thompson 1992, 
Paczynski 1992, Kouveliotou et al. 1998, 1999) create the highest energy cosmic rays (ultra- 
high energy cosmic rays, a.k.a. UHECR). I assume magnetars occur in all normal galaxies 
which form massive stars, with the UHECR ariving from outside our own Galaxy, except on 
the occasions (perhaps once per 10^ years) when a rapidly rotating magnetar is born in our 
own galaxy. This model has the virtue of having little difficulty in accelerating protons to 
energies in excess of 10^^ eV, with the sources being in all normal (star forming) galaxies and 
a luminosity density entirely acceptable from the point of view of the (very approximately 
known) rate of formation of magnetars in our own galaxy, thus offering an explanation of 
the puzzling air showers with energies above 10^^ *^ eV without having to introduce major 
extrapolations of known or suspected energetics. 

In §2, 1 summarize the data on UHECR, the loss processes that affect their transport 
through intergalactic space, the energetics of UHECR and aspects of the astronomy and 
physics of magnetized compact objects as known from studies of rotation powered pulsars 
that are relevant to the present investigation. §3 outlines the calculation of the particle 
injection spectrum from a newly born magnetar and the effects intergalactic transport have 
in altering the injected spectrum to the spectrum received at the Earth. In §4.1, I discuss 
the escape of the relativistic wind from a magnetar 's natal supernova, and in §4.2 I outline 
how the wind drives a relativistic blast wave containing a Magnetar Wind Nebula (MWN) 
into the surrounding interstellar medium. The accceleration mechanism of the UHE ions 
in the magnetar's wind is addressed in §4.3, and their escape from the expanding MWN 
nebula is outlined in §4.4. The effect the electromagnetic energy lost from a magnetar has 
in creating HI supershells in the interstellar medium is discussed in §5, with results that are 
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used to suggest that more than 90% of a newly born, rapidly rotating magnetar's 5 x 10^^ 
ergs of rotational energy gets lost as gravitational radiation. The typical strains of such 
gravity waves, from relatively nearby events that could contribute particles in the UHECR 
spectrum with energies above 10^° eV, are estimated in §6. The possibility of observing 
bursts of UHECR associated with the birth of individual magnetars gets attention in §7. I 
discuss the relation of this study to other work in §8. 1 draw my conclusions in §9. 

The most prominent difference between this model for UHECR and most others which 
appeal to sources more or less uniformly distributed throughout the Universe is that the 
underlying objects create an extremely fiat injection spectrum, with particles per unit energy 
range injected at a rate oc E~^[l + {E / EgY]'^ , with Eg reflecting the strength of gravity 
wave emission on the magnetars' early spin down; for the model developed in detail here, 
s = 1. As a result, the usual Greisen-Zatsepin-Kuzmin (GZK) cutoff disappears. This 
cutoff, apparent in models that assume power injection spectra oc E~^, s >2 (see Berezinsky 
et al. 2002 for a recent example of such "conservative" models), becomes a flattening of the 
observed E^J{E) spectrum between 10^^ and 10^° eV; the flatttening becomes a moderate 
replica of a GZK cutoff if the gravitational wave losses are strong enough. 



Recent observations of ultra- high energy cosmic rays (UHECR), manifested as atmo- 
spheric air showers with energy above 10^^'^ cV, suggest an origin in the metagalaxy (Bird et 
al. 1995, Yoshida et al. 1995, Takeda et al. 1999, Abu-Zayyad et al. 2002a,b); see Nagano 
and Watson (2000), Watson (2002) and Halzen and Hooper (2002) for recent reviews. This 
interpretation relies on three basic facts. 

1) The total cosmic ray energy spectrum hardens above Ea ~ 10^^'® eV, with isotropic 
intensity J{E) oc £'-3.2±o.05 i^gj^^ ii^^ to J oc £'-2-8±o.2 f^^^ ^ ^ hardening of the 

spectrum known as the cosmic ray "ankle." Sec Nagano and Watson (2000) for a thorough 
discussion of the strong evidence for the spectral hardening at the highest energies, including 
the uncertainties in the value of Ea. Such hardening of the spectrum indicates an origin for 
the high energy particles causing the showers at energies above Ea differing from the lower 
energy particle sources. Numerically, for E > Ea 



This flux yields the integral intensity above Ea, J{> Ea) ~ 3 x 10 cm ^-sec -"^-ster ^, 
corrsponding to particle number density n{> Ea) ~ 1.5 x lO^^'' cm~^ and energy density 
U{> Ea)^2x 10-^ eV cm-=^. 

2) Correlations of event directions on the sky with the supergalactic equator may exist 
at energies above 10^^'^ eV (Uchiori et al. 2000), although the results are also marginally 
consistent with isotropy, as are other analyses (Takeda et al. 1999 and references therein.) 



2. Background 
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An extragalactic origin remains consistent with either interpretation of the isotropy investi- 
gations. The few events with energy above 10^° eV say nothing useful about anisotropy. 

3) The shower data are consistent with hght nuclei or protons as the initiating events 
above At lower energies, in the cosmic ray "shin" , lying between the "knee" at 10^^ eV 
and the ankle, the data suggest a heavy nuclei composition (Bird et al. 1994). 

If light charged particles such as protons initiate the UHE showers, the Larmor radii of 
such particles in the interstellar magnetic field, = E/ZeB = 100(£;/1020 eV){10-^ Gauss/ZS) 
kpc, exceed the radial extent of the gaseous disk of our galaxy (for Z — 1), and greatly ex- 
ceed the vertical thickness of the disk - indeed, rL(100 EeV) greatly exceeds the thickness 
of the galactic synchrotron halo. Even ii Z = 26, as in the galactic pulsar source model 
of Blasi et al. (2000), the Larmor radius of 100 EeV particles still exceeds the thickness of 
the known magnetized disk (the synchrotron halo, with scale height ~ 1 kpc), implying a 
substantial reduction of the flux estimated by Blasi et al? The observed isotropy also argues 
against the UHECR being galactic or local group neutrons, which can survive for a distance 
~ 9(i?/10^^ eV) kpc from their sources; neutrons at 100 EeV certainly cannot arrive from 
the whole local supercluster. 

Extragalactic ultra-high energy charged particles cannot arrive from truly cosmological 
distances. As first noted by Greisen (1966) and by Zatsepin and Kuz'min (1966), interaction 
of relativistic charged particles with the cosmic microwave background (CMB) degrades their 
energy at a rapid rate through production of tt mesons when the CMB photons have sufficient 
energy in the rest frame of the cosmic ray. The attenuation length for protons due to pion 
losses (Protheroe and Johnson 1996, Stancv et al. 2000, Fodor and Katz 2000, Berezinsky 
et al. 2002) can be adequately represented as 
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Eg = 2.8 X 10^° eV, (2) 
Dgzk = cTgzk = 20 Mpc, Tgzk = 6.5 x 10^ years. (3) 



Expression (2) gives the largest loss rate above 4 x 10^^ eV. Below this energy, photopion 
losses are weaker than those due to pair creation and to cosmological expansion. For 
the purposes of the present study, pair creation and cosmic sxpansion can be adequately 



^From what is known of the top layers of a pulsar's atmosphere from X-ray observations, the composition 
almost certainly is not iron. Helium may be a more probable candidate, based on interpretaion of neutron 
stars' X-ray spectra [including line features in one case, Sanwal et al. 2002 and references therein, or perhaps 
O and Ne (Hailey and Mori 2002).] Lighter elements would have floated up to the X-ray photosphere (Zavlin 
et al. 1998, Pavlov and Zavlin 2000). Such high altitude ions (perhaps a few centimeters above the crust) 
are the reservoir for extraction from the star by the electric fields parallel to the almost vertical B fields, to 
be injected into the escaping wind, not the iron thought to dominate the crustal composition. Therefore, a 
charged particle model of the UHECR based on acceleration driven by neutron stars more plausibly should 
accelerate protons or a particles, as indicated by the shower structures at energies above Ea, which in turn 
requires consideration of an extragalactic model, to be consistent with the isotropy data. 
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represented as having constant attenuation lengths^, l± — 1000 Mpc, T± — 3.3 Gyr, and 
Ih — cTff — 4300 Mpc, Th — 14 Gyr. Therefore, the average loss rate is 



E Tioss{E) Th T± Tgzk 



Many hypotheses have appeared concerning the origin of these air showers. The most 
puzzhng aspect of the problem has been the observation of air showers at energies so large 
(in excess of 10^° eV) that charged particles coming directly from the accelerator must have 
started no more than 50 Mpc away, due to the radiation losses in the CMB. Within this 
local patch of the Universe, there is a lack of the most commonly hypothesized accelerators, 
disks and jets associated with giant black holes in the most active galactic nuclei, an idea 
with a long history going back to Lynden-Bell (1969). 

This scarcity has led to a plethora of explanations. Some are of the "bottom up" variety 
as in this paper, based on acceleration of normal charged particles to extreme energies, 
usually in objects not otherwise known to exhibit the ability to accelerate extremely high 
energy particles {e.g.^ Torres et al. 2002 and references therein). Others are of the "top- 
down" variety, based on the behavior of particles left over from the early Universe - see 
Sarkhar (2002) for a recent review of such ideas. The Z-burst model (Singh and Ma 2003 and 
references therein) offers a mixed alternative, with extremely high energy charged particle 
acceleration and ultra high energy neutrino production moved to high redshift (where hyper- 
active AGN are more common), while the high energy protons we observe are created locally 
through interactions between the UHE neutrinos and low energy neutrinos left over from the 
Big Bang (clustered in dark matter halos, if the background neutrinos have sufficient mass). 
The main point of the model proposed here is to argue that adequate accelerators may be 
common in the local Universe, since they may be in all normal galaxies with active star 
formation. In particular, the accelerators may be associated with magnetars, a relatively 
rare form of core-collapse supernova debris. 

The required galactic luminosity in UHECR gives a strong hint that the UHECR sources 
have something to do with massive stars and supernovae in ordinary galaxies. This luminosity 
is, using (1) and (4), 

i^r^ = — / dEE,^^ = 0.8 X lO^^n;,^ ergs/sec. (5) 

Here Ug — 0.02np2/Mpc^ is the space density of galaxies (Blanton et al. 2001), the majority of 
which actively form stars and are not in large clusters. The particles' intergalactic residence 
time, weighted by the energy in the observed spectrum in (5), is {Tioss) ~ Th/6.9 2.1 Gyr; 



^More accurately, l± = 500[1 + {E±/E)°-^] Mpc for E < E± = 4 x 10^^ eV, but the slow increase in 
path length with decreasing energy has no effect on the observable aspects of the model outlined here, at its 
present level of development. 
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when weighted by particle number, the mean residence time is ~ Th/3.8 = 3.7 Gyr. This 
luminosity per galaxy coresponds to the metagalactic luminosity density 5 x 10"^^ ergs/Mpc^- 
year. These energetics are consistent with a supernova related source, as is the case with 
normal cosmic rays. The high energy per particle forms the distinctive peculiarity of the 
UHECR. 

The basic model for the accelerator proposed here is a variant of that advanced by 
Blasi et al. (2000), who assumed the relativistic winds from pulsars in our own galaxy 
could be the acceleration site, granted that these stars begin their lives as almost maximally 
rapid rotators. Relativistic magnetic rotators have magnetospheric voltage drops across the 
magnetic field with magnitude 

= ^ = *mag- = 3 X 10-—^ f T7#3t)' Volts, (6) 
c lO'^"^ cgs \ 10^ s ^ / 

which is in excess of the energy / charge exhibited in the UHECR, if the underlying star has 
large angular velocity {VL > 10^ s~^) and large dipole moment {fi > 10^^ cgs, corresponding 
to a surface dipole field above 10^^ Gauss)"^. Here ^mag = R\B{Rl) is the magnetic flux in 
the "open" field that extends beyond the light cylinder, located at r = = c/fl = 30/1^4 
km that connects the magnctosphcrc to the outside world, and Q = 10^1^4 s~^. 

Acceleration in the wind avoids the catastrophic radiation losses associated with accel- 
eration inside a rapidly rotating neutron star's magnetosphere. Gharged particles moving 
on paths with radius of curvature pc have radiative power P = {2/3){q^c/ pl){E /mc^)'^. Ra- 
diation losses would therefore prevent protons from accelerating to energies above 
Erad ~ lO^^-^{A^/Zy/'^Q^^^\r/RLy/%Pc/r)^/^, eV, where a particle's mass is ArUp and its 
charge is Ze. Within and near the magnetosphere, the strong magnetic field guides particles 
along B. Therefore, even if strong parallel electric fields can form, the radius of curvature 
of a particle's orbit is Pc ~ ~ -^l and all the energy expended in accelerating the putative 
cosmic ray gets lost as gamma rays (and pairs.) 

At radii r » Rl, the physical situation differs. I shall argue that at large radii, acceler- 
ation can and does occur in an angular sector around a magnetar's rotational equator, where 
the radii of curvature of highly accelerated ions' orbits do greatly exceed the local radius r 
at each point on their orbits (see §§4.3 and 4.4), so that radiative losses are negligible. 

Relativistic magnetic rotators, appearing as rotation powered neutron stars and as black 
holes with magnetized disks losing rotational energy in jets, generally create non-radiative 



''indeed, a neutron star which had magnetic energy in the poloidal field comparable to the gravitational 
binding energy of the star (Sdipoio ~ 10^* Gauss) would be able to accelerate charged particles to ^ 10^'"' eV, 
if it rotated close to breakup > 10^ sec~^). Of course, such an object would look nothing like a sphere 
and would lose most of its rotational energy in gravitational waves (see below), should it be so fortunate 
to survive the host of instabilities plaguing its existence. Magnetized fallback disks {e.g., Macfadyen and 
Woosley 1999) might be functional equivalents of magnetically supported, rapidly rotating neutron stars; 
see §8. 
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relativistic outflows. These outflows - relativistic "winds" - include wound up magnetic flelds 
and outflowing plasma, which support the "MHD" electric field Ej^ = — (l/c)v x B. The 
magnetic field has magnitude B ^ B{Ri){RL/r), an estimate appropriate both to quasi- 
spherically expanding winds and to (current-carrying) jets. At radii large compared to R^ , 
the cross field voltage in the wind is (v/c ~ 1) 

*«,.nd ^ rE^ir) = rB{r) = RlB{Rl) = ^ ^ ^ (7) 

The last two forms of (7) are appropriate for underlying rotators which have magnetic dipole 
moments, such as neutron stars, while expressing ^wind in terms of the open magnetic flux 
^ allows application of the same physical considerations to outflows from the magnetized 
disks thought to orbit black holes. 

The estimate given by equation (7) assumes the relativistic wind does not suffer strong 
radiative dissipation, which is true in the outflows from rotation powered pulsars, and is 
also the case in electromagnetic models of jets from underluminous galactic nuclei in radio 

galaxies, galactic microquasars and gamma ray bursts (GRBs). Then the voltage generated 
by the magnetized rotator is available to accelerate particles at radii r ^ c/Vt - for protons, 
particles of energy greater than lO^'''^ eV could be accelerated at distances r > 0.3(r/pc)^ 
pc, a dimension easily achievable in the winds from magnetized rotators. If the acceleration 
process approximates that of a linear accelerator, rather than of a synchrotron, as may be 
the case if "surf-riding" acceleration (e.g., Tajima and Dawson 1979, Chen et al. 2002, §4.3) 
applies, then r/ 1 and radiation losses are even less of a constraint on the model (see 
§4.4). 

The new aspect of the idea pursued here is to apply the same physics to magnetars in 
all galaxies. These objects have no difficulty (in principle) with accelerating light nuclei to 
the required energies, and by hypothesis, being located in all galaxies, have no difficulty in 
supplying particles with an isotropic distribution of source directions. The initial spin rates of 
magnetars are unconstrained, except from theoretical considerations of the stars' dynamical 
stability, which limits the rotational energy to be less than 25% of their gravitational energy 
[Vt < 1.2 X 10^ s~^ for a uniform sphere). 

Therefore, one can assume that the stars are born with sub-millisecond initial rota- 
tion periods. The limited evidence on the initial spins of rotation powered galactic pulsars 
suggests some begin their lives rotating at less than their maximal rate, z/ electromagnetic 
torques dominate their spindown from birth (e.g., Murray et al. 2002, Kaspi and Helfand 
2002). As I will show, other constraints suggest that at least for magnetars, gravitational 
wave torques may govern the initial spindown, which alleviates concerns over whether all 
such objects might be born as slow rotators - in any case, only 5-10% of the magnetars must 
be born as rapid rotators, for the UHECR source model proposed here to be viable. 

Magnetars are of particular interest as an acceleration site since a) they exist (Kouve- 
hotou et al. 1998, 1999; see the summary in Baring and Harding's 2001 Table 1), and b) 
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if they rotate fast enough sometime in their hves, they easily can have voltage drops well 
in excess of the energy per unit charge observed so far in UHECR (expression 6). ^4 = 1 
corresponds to a rotation period of 0.63 msec, about half the ccntriftigal breakup angular 
velocity for a IAMq neutron star with 10 km radius, while ^^33 = 1 corresponds to a polar 
surface dipole field of 2 x 10^^ fi^sRiQ Gauss, with Riq = km; these values corrspond 

to a charateristic voltage of 3 x 10^^ V. Clearly, rapidly rotating, magnetized neutron stars 
can provide the energy/particle seen in UHECR, in principle, so long as the wind carrying 
the voltage survives to distances far from the star. 

The energetics of a magnetar acceleration model are simple. Suppose the particles each 
gain the energy 

E{Q) = g77$^ag = = 3 X lO''^ ZrjiQlfi^s eV, 771 = rj/O.l. (8) 

Here r] is the fraction of the open field line voltage experienced by each particle on its way from 
the star to the outside world. If the magnetosphere is relativistic, the usual electromagnetic 
rate of rotational energy loss 



applies^. 



I adopt the view that electric conduction currents are a major part of the support 
for the electromagnetic fields which exert torques on the rotating objects. I also adopt 
the assumption that charged particle beams extracted from the rotating object by electric 
fields compose these currents; see Arons (2002) for a recent discussion of these currents' 
significance, and a discussion (briefiy recapitulated below) for the evidence that such beams 
actually exist. 

If conduction currents are the sole source of the electromagnetic fields, and if a charged 
particle beam carries all of the current, Er — /$mag — qcN^xaag-i and 



N 



dipole 

27rc 



\q\c 



(27rAap), (10) 

(11) 



mag 



One recognizes Acap'n'i^R^/cjRl to be the area of a magnetic polar cap, if the closed field 
lines extend to the light cylinder, and the expression in square brackets in (10) to be the 
Goldreich- Julian charge density pcj at the stellar surface, the number of elementary charges 
per unit volume such that parallel electric fields are shorted out in the magnetosphere. 



^Expression (9) is the standard result for vacuum dipole radiation, after averaging the geometric factor 
(2/3) sin^ i over the sphere; i is the angle between the dipole moment and the rotation axis. All electromag- 
netic spindown theories are thought to lead to comparable rates of energy loss, with Eem = c^mae- 
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The electric current 

Igj = 2A,apCpGJ = = C$„ag (12) 

is a fundamental item in the theory of the spindown of relativistic magnetized rotators. Since 
Eem — loj^mag, rotators which have such currents linking the star to the outside world can 
spin down at a rate comparable to that expected from vacuum (magnetic dipole) radiation 
losses. Both contributions to spindown are expected to be comparable, in oblique rotators - 
there is substantial evidence, from the lack of dependence of radio pulsars' spindown torques 
on the angle between the rotation and dipole axes (Lyne and Manchester 1988) that this 
assumption is true. 

Such currents, in stars with dipolar magnetic fields, include a primary charged particle 
beam emanating from the polar cap - electrons, when i = Z(/x, ft) < 90°. That current 
must be balanced, at least on average, by an electric return current. The most traditional 
ideas assume this return current is an outflowing beam of charges of the opposite sign - ions, 
if i < 90° - usually assumed to occur on the ring of field lines around the polar cap which 
form the boundary layer separating the closed and open field line zones of the magnetosphere 
(Goldreich and Julian 1969); in force free models (Michel 1975, Contopoulos et al. 1999), the 
return current appears as a current sheet. Such "auroral" field lines map to the rotational 
equator, in oblique as well as aligned rotators, where the return current forms a (corrugated) 
current sheet, with spatial oscillations imposed by the rotation of the underlying non-aligned 
dipole magnetic field. The model developed here for UHECR assumes the ion injection rate 
per magnetar to be that of this electrodynamically fundamental return current, based on an 
analogy to what has been learned from Pulsar Wind Nebulae. 

Eem = Icj^open does not mean that such an object automatically is a maximal particle 
accelerator, with all the energy going into the energy of the particles in the electric currents. 
For that to be true, the particles in the currents would have to transport all of the rotational 
energy lost, and each would have have to reach the maximum energy per particle g$mag - 
other particle constituents would have to be largely absent, and Poynting fluxes would have 
to be unimportant. Pulsar theory, and related theories of magnetized disks, begins with the 
opposite assumption - that a large scale Poynting flux carries almost all of the outlowing 
energy, with the current carriers gaining very little of the potential energy available until 
dissipation occurs in the region where the outflow from the star terminates (Goldreich and 
Julian 1969). Even then, the energy dissipated might go into bulk expansion of the sur- 
roundings, rather than into acceleration of a few high energy particles. This approximation 
- represented in its most extreme form by force free models of the outflow {e.g., Contopoulos 
et al. 1999) - gets its emprirical support from the lack of intense photon emission from 
many known relativistic outflows. Therefore, one cannot blithely assume the existence of an 
energy /particle given by (8) with ~ 1, even when the ion loss rate is given by (11). 

In addition, the study of pair creation in the magnetosphere has led to the belief that a 
dense electron-positron plasma provides the dominant constituent by number in the outflow. 
The feeding of the X-ray source in the Crab Nebula, in the X-ray nebula around the Vela 
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pulsar (Helfand et al. 2001) and in the nebula G320.4 around PSR 1509-58 (Gaensler et al. 
2002) strongly supports the existence of such a dense plasma outflow - if the flow leaving 
the neutron star has the same composition on all field lines, the particle input to the Crab 
(A^totai ^ 10^^ pairs/second) suggests i]pairs < 10^^ - the pairs carry much of the energy flow 
at large distances from the star, but being a dense plasma, they have energy per particle 
well less than (8). 

However, study of the structure of the shock wave terminating the equatorial relativistic 
wind in the Crab {e.g. Hester et al. 1995), in G320.4 (Gaensler et al. 2002) and perhaps in 
the Vela nebula (Helfand et al. 2001), where the energy in the pulsar wind gets transfcrcd to 
the plasma forming the observed synchrotron nebulae, suggests that an interestingly large 
fraction of the total voltage does get apphed to a Goldreich-Julian flux of heavy ions on the 
equatorial flow lines from rotation powered pulsars (Arons 2002). These are probably hght 
nuclei such as protons or alpha particles (Gallant and Arons 1994). Detailed modeling of 
this interaction zone (Gallant and Arons 1994, Spitkovsky and Arons 2000 and submitted 
to ApJ) suggests that the total ion injection rate is indeed the Goldreich-Julian rate (11) 
and that these these particles, flowing out in a latitudinal sector around the stars' rotational 
equators, have experienced ~ 10% of the full potential drop (i.e., rj ~ 0.1), and thus carry 
an energy loss from the rotator Eions = lajV^mag = V^R ~ 0.1£^r of each neutron star's 
rotational energy loss in an equatorial sector filling ~ 1/5 of the sky around each of these 
pulsars, exactly the region where the models hypothesize the fiow of the return current in 
an equatorial current sheet. Outside this equatorial sector, the wind is more likely to be 
Poynting fiux (AC and DC) dominated. 

The specifics of the ion accelerator are not known - most recent suggestions for the 
acceleration site have focused on non-ideal processes in the equatorial wind, well outside 
the light cylinder. The ions enter the wind, having been drawn up from the star along the 
boundary layer separating the closed and open field lines, leave the star in the plane of 
the star's rotational equator, and gain their high energy per particle at some large distance 
from the star, perhaps at the interesting radius where the shortage of charge carriers causes 
dissipation of the corrugations in the equatorial current sheet (Michel 1971, Coroniti 1990, 
Michel 1994), or transform to electromagnetic waves, either Alfvenic (Bellan 1999, 2001) 
or vacuum-hke (Melatos and Melrose 1996; Melatos 1998), whose pomderomotive force can 
accelerate the particles, through surf riding in the EM fields, to energies comparable to 
Q^mag- I discuss surf riding acceleration mechanism further in §4.3, but a detailed study of 
that physics is outside the scope of the present investigation. 

In what follows, I assume that such ion acceleration also occurs in magnetars that were 
rapidly rotating at birth, and explore the consequences of this assumption for an extragalactic 
origin of UHECRs from normal galaxies. I assume all galaxies which have core collapse 
supernovae form magnetars with a birthrate comparable to that in our own galaxy, 
where 10~^ > I'm > 10~^ years"^ (Gaensler et al. 2001), and that some are born with 
high angular velocities, fli ~ 10^ s~^. During their spindown, they emit heavy ions (taken 
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to be protons in what follows) with an ion loss rate per magnetar given by (11), with the 
energy/particle in the emitted equatorial "beam" declining as the rotator ages, as specified 
by (8), with 77 «i 0.1 in analogy to what has been inferred for pulsars. Therefore, the highest 
energy particles come from new born objects. The particles cannot be contained within 
normal galaxies, therefore the flux at the earth reflects the contributions from all galaxies 
within the energy dependent GZK volume. That volume expands to encompass all galaxies 
in the Hubble volume, as the particle energy decreases, but below the ankle energy this 
extragalactic source becomes swamped by the softer particle spectrum thought to be due to 
acceleration in our own galaxy. 

A simple estimate shows that the idea is a viable candidate for the origin of the UHE 
spectrum. The star emits ions at the rate (11). Suppose electromagnetic torques are the sole 
means of spinning down the neutron star. The strong Q dependence says that the star emits 
most of the particles during the initial loss of rotational energy, when Q — 0,^^ — 10^Q4. The 
initial electromagnetic spindown time is 

Tem = -rr- = ^-TPyl = 5 -TTyl mmutes. (13) 

Therefore, each magnetar injects 

N, ^ Noj^TEMm = = 4 X 10^2^ (14) 

ions into the metagalaxy, each with the energy (8). The existing observations of a few events 
at energies > 10^° eV require the birth of some magnetars with rjiQli^ss > 0.1/ Z. 

With a magnetar birth rate per galaxy i/^ = lO'^^Urn^ yr~^ and a galaxy density Ug = 
0.02ng2 Mpc~^, the number of particles per unit volume from magnetars in all the galaxies 
is 'n'uiiE,magnetars = NingymTioss{Ea) , where 7]oss(-E^a) IS the residence time of a particle at 
the ankle initially injected at higher energy; very roughly, this is the hfetime T±{Ea) ~ 10 
Gyr. Then 



J^'''\> Ea) - ^NiUmUgT^iEa) ^7x 10"^^ (^) cm-^ - sec-^ - ster-\ 
comfortably above the observed J(> Ea) = 3 x 10~^® cm~^ — sec~^ — ster~^. 



(15) 



The required luminosity in UHE particles does not unduly tax the energetics of rapidly 
rotating magnetars in galaxies, as can be seen by writing L^^^ — Smi'mE^ . With Lg^^ 
from (5) and E'j^^ = {l/2)mf = 5 x 10^2j^^^2 ^^g^^ ^^^^ 

require a fraction e„j = 

0.005/z/m4^g2-^45^4 of the initial rotational energy to go into UHE particle acceleration. I 
show in §5, using the results of the more detailed theory in §3, that the birth rate of fast 
magnetars, those with initial voltages large enough to accelerate particles with energy greater 
than is 1/^"**'' = (0.05 — O-Vju^. Then the coefficient in changes from 0.005 to .05-0.1, 
still not a taxing demand on the model's energetics. Put slightly differently, the luminosity 
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density in rapidly rotating magnetars deduced from the fit of the model to the observations 
described in §3 is - v^i^'^^ngE^^ = 10^^(z/^"'*VlO"^ yr~^)h5^l ergs/Mpc^ - year. Given 
that 10% of the accelerating voltage gets applied to each accelerated particle, and that the 
injection rate per object is given by (11), the model has the required luminosity density 
~ 10^^ ergs/Mpc^-year in UHECRs. 

As I shall show, more careful treatment of the spectrum, plus approximate accounting 
for the losses sustained as the magnetized wind bursts out of the supernova that gives rise to 
the neutron star, improves the correspondence between the model and the data, and makes 
some interesting predictions for future UHECR observations. 



3. The Spectrum of UHE Cosmic Rays from Bare Metagalactic Magneteirs 

3.1. The Injection Spectrum from an Isolated Magnetar 

During a time dt, the magnetar spins down from angular velocity Q{t) to fl{t + dt) — 
il{t) + Cldt. According to (8) and (11), the energy/particle E{t) — miC^'y{t) and the instan- 
taneous injection rate TVj each decay oc Jl^. Therefore, the number of particles accelerated 
as the star spins down in time dt with energy between mc^7 and mc^(7 — ^7) is 

-dNi = Ni[n{t)]dt = Ni^^dj, (16) 

dil 07 

with 

_ A _ Eji _ Eem + Egrgy , , 

dt ~ ~ m~ m ' ^ ^ 

Here / is the moment of inertia, and (9) expresses the electromagnetic energy losses. 

The gravitational wave losses may be substantial. If the star has a static equatorial 
ellipticity e, gravitational waves extract energy at the rate (Ostriker and Gunn 1969) 

32 GPe^Vt^ 

-Egrav = -T -5 = 1-8 X lO^'ll^lel ergs/sec, (18) 

c 

where = //lO*^^ g-cm^ and 62 = e/10~^. Note that €2 ~ 1 is a substantial but not 
ridiculous ellipticity; interior magnetic fields can create ellipticity (Usov 1992, Bonnazola 
and Gourgoulhon 1996) with magnitude (Ostriker and Gunn 1969) 

_ ^5 R ^gp2 _ /rj2 \\ ^ ~ (■^toroidal) , /.qx 

2AGM^^ ^^toroidai))- (4 X 1016 Gauss)2 ' ^^^^ 

the numerical value assumes M = IAMq. The initial spindown time due to gravitational 
radiation is 

yn^i 28.5 , , , 

tgw = —r = 7 — 2m seconds. (20) 
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Interior fields an order of magnitude stronger than the surface field are by no means 
incredible. Formation of interior fields as large as 10^^'^ Gauss in magnetars has been pro- 
posed, as a consequence of dynamo activity in the first few seconds of the proto-neutron 
star's life (Duncan and Thompson 1992, Wheeler et al. 2000.) - these stars are thought 
to have non-dipolar surface fields much stronger than the already large dipole field. Once 
formed, such an ellipticity would most likely persist throughout the neutron star's spindown, 
certainly through the first day or so of the spindown germane to the model discussed here. 

Other sources of ellipticity or its equivalent may be dynamical instability of the star, 
if the rotational energy exceeds 25% of the newly formed star's gravitational energy, and 
secular instabilitcs of F- and R-modes, driven by the gravitational radiation itself (e.g., Lai 
2001, Ushomirsky 2001). In the interest of simplicity, I use (19) as motivation to choose 
£2 ~ 1 as a parameter throughout the subsequent discussion, and neglect other versions of 
gravitational wave spindown. In the appendix, I briefiy summarize the possible effects of 
gravity waves emitted because of the R mode instability. 

I rewrite the spin down rate as 



-n 



9 /c3 



a; 



(21) 



with Vtg the angular velocity at which gravity wave and electromagnetic losses are equal. 



5 cV" 
72 GPe'^ 



1/2 



= 3 X lO^i^ s-\ 

J 45^2 



(22) 



Expressions (8), (11), (16) and (21) then yield the particle spectrum accelerated by a mag- 
netar during its spindown to be 

-1 



dNi dNi 



d'j 



dt 



dt 

In 



mc 



,dQ 9 c'l 



dE 



where 



Eg = niiC^-fg 



Zrjefj, 



4 Ze/ij 
5 Z-qeji^ 



1 + 



72 GPe^ 



3 X 10 



20 



2 

t2 2 
-'45^2 



eV. 



(23) 



(24) 



The star promptly emits the energy that goes into the wind. The time to spin down to 
a voltage such that a particle of energy E can be accelerated is 



tspin{E) _ I ^ ] _ [ iiy ] 



^-^-In 
E Ei 




E 1 + 





Q2" 


















Eg J 





where 



2-'"g 
4M'f^| 



81 GPce^ 
~5 



0.9 



j3 2 
- '45^2 

f^33 



hours 



(25) 



(26) 



(27) 
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is the electromagnetic spin-down time, for a star spinning at angular frequency Jig. (13) 
and (20) give the initial spin-down times for pure electromagnetic and pure gravity wave 
losses. Expression (25) yields the spindown history when gravity wave losses dominate, 
Q(t) = Qi[l + {2t / Tgw)]^^^^ J if the limit Qg <^ Q < Qi exists, and yields pure electromagnetic 
spindown, Q{t) = + (t/rgM)]"^''^, in the limit Q -C ^g- For t < 2Tg, (26) simplifies to 



tspin{E) — 2tgw 



- 1 

n 



1.8^^^ (^^^-^ ) hours, < E < (28) 



Here E^ax = Zerj^^aagi^i) = 3.3 x 10"^^ Zrjifi^s^l eV. For t > 2Tg, E < Eg, electromagnetic 
losses alone accurately describe the spindown. In this limit, 

(02 \ 7 T 1 r|20.5 „\r 

g - 1 j 1-8 ^2^' E hours, E < Eg. (29) 

The general spin history can be written in the simplified form, valid for t ^ tqw, 

^ = ^g(^Y, (30) 



t 

with p = 1/4 for tqw <^ t < 2Tg and p = 1/2 for t > 2Tg. 

Expression (23) is the spectrum that would be observed from a single nearby event 
by an experiment which accumulates events of energy E and higher, over the time (26). 
An experiment which resolves a single nearby magnetar birth and spindown would see an 
instantaneous, approximately monoenergetic spectrum, with the energy/particle declining 
in proportion to Q'^(t). 

To the magnctar's host galaxy and to the metagalaxy, the energy input to the relativistic 
wind appears as an impulsive burst; for the energies of all the observed particles {E^ax ^ 
= 6 X 10^^ eV), tspin{E) < tspin{Ea) = 42Z771/45///33 hours. Note that this is the 
time taken to put energy into the wind, corresponding to voltages high enough to accelerate 
the observed ultra-high energy cosmic rays; the actual time of acceleration is later. How 
much later depends on the time it takes for a region of the wind to flow to the radius where 
acceleration occurs. In §4.3, I suggest this flow time probably is a few hours, for typical 
parameters. 



3.2. Intergalactic Particle Transport and the Observed Spectrum from Bare 

Magentars 

3.2.1. Scatter Free Transport 

The magnetars forming in galaxies emit UHE particles which cannot be contained within 
the galaxies. Therefore, normal galaxies inject particles into intergalactic space at the average 
rate per unit volume Ugi'^dNi/d'j. Assuming scattering in intergalactic magnetic flelds has 
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negligible effect on particle transport, the intergalactic spectrum A'"(7) may be determined 
from 

d dNi 

— (^AT) = WgeomT^mrig-^. (31) 

Here 7 = —■j/Tioss{l) and Wgeom = 0.5 represents the geometric fact that only half the 
magnetars are expected to have the magnetic geometry appropriate to injecting ions from 
the stars atmospheres into the wind in the rotational equator (§4.4). 

Integrating equation (31) with (4) and converting from A^(7) to J{E) through rriiC^J = 
cN/An yields 

J(E) = ^_iV(7) = W,,om-^ ^""^ ngUmTioss(E) f Ema. ^ + . . . 

^ ^ inrriiC^ ^'^ ^^"^IGnZefx E y E 1 + ' ^ ^ 

The spectral form has the 1/E shape of the injection spectrum below 10 EeV, steepens 
to become proportional to E~^ in the observed UHE range 10 EeV < E < 100 EeV, then 
flattens to return to a shape approaching the E'^ spectrum above 300 EeV (oc E^^ if 
gravitational wave losses during spin down truncate the number of particles produced at 
energies well below E^ax, i-^-, when Eg <^ E < E^ax), where GZK losses saturate. In 
the strong GR case, gravity wave emission shortens the time at which each magnetar can 
contribute particles at a given energy, which reduces the flux even in the better studied 
region £■ < 5 x 10^^ eV, for parameters other than ellipticity fixed. 

These ultra high energy particles have little trouble entering our galaxy^. Therefore, 
expression (32) can be compared directly to the observations. Multiplying (32) by E^ and 
overlaying the results on the data yields the result shown in Figure 1. Given the substantial 
errors, fitting by eye suffices, with the conclusion that this elementary theory provides a 
satisfactory fit to the data for the values of Kq — Wgeom^mliT-gTH / pi- shown in Table 1. 



'^The ion Larmor radius = E/ZeB = 2.3{E/10'^^-^ eV)(3 ^Gauss/ZBo exp(z/2/i) kpc, where z is the 
vertical distance from the galactic midplane and Bq is the midplane magnetic field, is larger than the vertical 
extent of the synchrotron halo of our galaxy, whose scale height is /i ~ 1.8 kpc (Beuermann et al. 1985) 
implying easy entry even at the low end of the UHE spectrum. An extended, quasi-spherical magnetized halo, 
if it exists, probably has a magnetic field at least an order of magnitude weaker (with Larmor radius ^ 20 
kpc), an assumption supported by radio synchrotron emission observed in other disk galaxies (c.g, Hummel 
et al. 1991, Lisenfeld and Yolk 2000). Thus even at the lowest observed UHECR energies, individual 
particles can reach our galaxy's midplane following almost rectilinear orbits, in times not much exceding 
the light transit time - at the lowest energies, the guiding center drift velocity might drop a bit below c, to 
~ c{rL/Rhaio)j if an extended halo with radius larger than 20 kpc exists. Since energy losses are negligible 
on this scale (and they would be negligible, even if the particles had to diffuse in the galaxy's magnetic field), 
the particle density and spectrum inside our galaxy comes to equilibrium with the intergalactic spectrum 
at the same density and spectrum as is in the metagalaxy (Ginzburg and Syrovatskii 1964, Longair 1994). 
However, magnetic deflection in halo magnetic fields probably does corrupt the accuracy with which one 
can infer a possible source's location from the look back direction along an air shower (itself known to ^ 1° 
accuracy) for energies below 5 x 10^® eV. 
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Fig. 1. — Comparison of bare magnetar theory to UHECR data. The experimental data are 
from Nagano and Watson's (2000) summary except for the data from the Hi-Res (monocular) 
results, which are taken from Abu-Zayyad et al. (2002b). The three theoretical curves are 
for differing assumed equatorial ellipticities of the neutron stars: e = ("no GR"); e = 10^^ 
("moderate GR"); e = 10~^ ("strong GR"). The theoretical curves were fit to the data by 
eye, by requiring the model flux at 3 x 10^^ eV to provide the whole intensity at this energy. 
That fit yields values for the average source strength Kq = Wgeomi^mlngTH / which are 
fisted for the three cases shown in Table 1. All the curves use the same upper cutoff energy 
Emax — 3.3 X 10^^ eV, corresponding to ZriiiJL^2,^\ = 1 (a polar dipole field strength of 2 x 10^^ 
Gauss and an initial rotation period of 0.6 msec.) Lowering the upper cutoff to the minimum 
acceptable value of 3 x 10^° eV yields qualitatively unacceptable fits to the observations - 
there is too much curvature in E^J{E) below 6 x 10^^ eV. Including a "galactic" component 
J gaiactic{E) = 7 X 10^'^(£'/30 EeV)^'^'^ ev^/meter^-sec-ster constructed to represent the 
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At the relatively low energies now observed, the presence or absence of gravitational wave 
losses makes little difference unless the gravity wave losses are large. Below 4 EeV, this ex- 
tragalactic particle spectrum hides below the flux of (probably) galactic particles at energies 
below the ankle. 

The spectrum in Figure 1 represents an average over many particle bursts, each lasting 
a couple of days (for E > Eg). These bursts should be distinguishable as such, with sufficient 
collecting area. One can readily show, using the time dependent form of (31) with its right 
hand side set to zero and with (23) as the initial condition, that the fiuence JF of a single 
burst has the same spectral form as is shown in Figure 1, but with the upper cutoff energy 
now determined by the losses in the CMB degrading the maximum injection energy to the 
value determined by the source's distance. 

Because the injection spectrum is extremely fiat (oc E~^ , if gravitational wave losses 
are negligible), the losses in the CMB flatten E^J but do not produce the GZK pileup 
characteristic of steeper injection spectra {e.g., Berezinsky et al. 2002.) If the claimed 
observation of such pile up by the Hi-Res experiment (Abu-Zayyad et al. 2002 a,b) is 
confirmed and the discrepancy with the AGASA results (Hayashida et al. 1994, Takeda et 
al. 1999; see http://www-akeno.icrr.u-tokyo.ac.jp/AGASA/ for results through Fall 2001) 
is resolved, the existence of such a feature would, within the context of the model desrcibed 
here, point to very rapid early spindown due to gravitational wave losses, which has the 
effect of steepening the highest energy part of the injection spectrum - see the "strong GR" 
curve in Figure 1. 





e 


E, (eV) 






no GR 


< 0.003 


oo 


0.03 


0.02 


moderate GR 


0.01 


3 X 1020 


0.06 


0.04 


strong GR 


0.1 


3 X 10^^ 


1.25 


0.81 



Table 1: Source strength Kq = Wg^omiJ^JlO-^ yr'^) (7^^/0.02 Mpc-^) {h^ / fis3)(TH/U Gyr) 
for the theoretical curves shown in Figure 1. /ST^^^^^ is the source strength inferred if one 
extrapolates the spectrum found at energies well below Ea into the UHE regime, then su- 
perposes the model. 



3.3. Predicted Spectrum at Higher Energies 

Unless current experiments are so fortunate as to have observed the upper cutoff, the 
spectrum of high energy particles should continue into the higher energy region where current 
detectors lack sensitivity. The specific model used here suggests the spectrum E^J{E) should 
resume rising with E above the energy Eg = 2.8 x lO^'' eV where the GZK loss rate becomes 
approximately energy independent (unless the gravitational wave losses are large, as shown 
in Figure 1.) 
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This prediction is simplified, in its assumption that all magnetars have exactly the same 
starting voltage (lO^^'^ Volts) and single particle efficiency rj ~ 0.1. The inferred value of 
Kq in Table 1 being small compared to unity in the weak and moderate GR cases has its 
most natural interpretation in the initial voltages (initial rotation periods, if all objects have 
the same dipole moment) having a distribution of values, with only a small fraction having 
voltages such that particles of energy greater than Ea can be accelerated. Since the results 
of §4.4 yield Wgeom = 0.5, Table 1 suggests that the birth rate of magnetars with voltages 
high enough to contribute to the spectrum in Figure 1 is only 6% of the overall magnetar 
birth rate, in the "no GR" case (a "fast" magnetar formation rate of once every 170,000 
years/galaxy); for "moderate GR", the required rate is 10% of the overall magnetar birth 
rate (once every 10^ years/galaxy). 

A distribution of initial voltages can have potentially observable effects on the predicted 
UHECR spectrum. To illustrate, suppose the distribution of birth rates as a function of 
starting voltage is a power law. 



with ^rnvn < < ^i,max, ^i,max = 3 X 102^33(^^^/10^ SeC-^)^ and ^rnvn = lO^'' VoltS, 



the potential drop characteristic of the known magnetars - as with pulsars, initial rotation 
periods comparable to the observed periods are possible, so that initial rotation periods as 
long as 10 seconds must be considered. Expression (32), with dvrn/d^i from (33) replacing 
Vm-, now gives E^dJ{E;^i)/d^i, with 7max replaced by 7^ = Zerj^i/mc^. Integrating over 
all initial voltages has the effect, when s > 1, of replacing the logarithmic cutoff in (32) 
by {E/Ei^max)^^'^^^^ S{E/ Ei^rnax, s), whcrc 5" is a slowly varying function which produces a 
more complicated logarithmic cutoff of the spectrum at Ei^rnax = Zer]i<^i^rnax- If s < 1, there 
is no substantial change from the single initial voltage model - when the spectrum of birth 
potentials is flat (s < 1), essentially all the stars would have starting voltage equal to ^i,max- 

Fitting the resulting modified spectrum to the data (again by eye) yields the required 
reduction of the normalization if s f=i 1.2, in the "no GR" case. The extra E~^ '^ now 
introduced in E^J{E) at energies below 10^° eV is undetectable in the current experiments, 
but with increased data above 5 x 10^^ eV the shape of the spectrum will provide significant 
constraints on such refinements of the model. Above 10^° eV, the decreased number of 
contributing sources can substantially reduce the high peaks at high energy shown in Figure 
1, by as much as a factor of 3. However, the basic conclusion, that the model predicts 
something to see, roughly at the level one would guess by extrapolating the observed E^J{E) 
above the formal GZK cutoff, remains robust. With s = 1.2, 5% of all the magnetars have 
Zerj^i above Ea, consistent with the simple inference of the fast magnetar birthrate described 
above. 




(33) 



Also, the highest energy part of the spectrum probes acceleration in the wind at the 
earhest times. As discussed in §4.1, breakout of the wind from a magnetar's natal supernova 
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dissipates a fraction l—Wuowout ~ 0.3 of the initial energy injected, with voltages in the freely 
expanding wind available for particle acceleration only after spindown from to ^jW^f^^touf 
suggesting a more realistic upper cutoff for the spectrum to be at EmaP — Emax^u^out ~ 
8 X 10^%i^33nl{Wuowout/0.7y/^ eV. 



Figure 1 illustrates another robust conclusion, that the highest energy part of the spec- 
trum is sensitive to possible spindown losses due to gravity waves. The shortening of the 
time in which particles can be accelerated to energies above Eg suggests the interesting pos- 
sibility that observations of UHECR at energies above 10^° eV can point the way to sources 
of gravitational radiation, and vice versa (see §6). 

While the agreement of this simple model with the observations is encouraging, one 
must consider the possible modifications due to losses. Blowout from the source (in this case, 
the supernova initially containing the relativistic wind), radiative and adiabatic expansion 
losses are the most prominent pitfalls for any compact object model for high energy particles 
observed in any part of the cosmic ray spectrum. Likewise, modification of the spectrum due 
to possible diffusive transport in the intergalactic medium must be considered. I consider 
scattering in intergalactic magnetic fields first. 



3.4. Effect of Pcirticle Scattering in an Integalactic Magnetic Field 

There is no positive information pointing to the existence of a general intergalactic 
magnetic field. Various hmits suggest the magitude of such a field does not exceed ~ 10~^ 
Gauss {e.g., Barrow et al. 1997), although a number of authors have argued for larger fields 
injected from galaxies and spread around by ill- understood processes {e.g., Farrar and Piran 
2000, Kronberg et al. 2001); see also Eilek and Owen (2002). If such a field exists, and if 
it causes high energy particles to flow diffusively through intergalactic space (Adams et al. 
1997), significant alterations of the UHE spectrum received from sources in galaxies would 
occur. 

Consider what happens in a magnetic field strong enough, and with sufficient magnetic 
fluctuations, to cause particles to have a scattering mean free path small compared to the 
GZK length at energy E. Assume the intergalactic field and its fluctuations flll all space (no 
intermittency). The mean free path for a particle of energy E in the fluctuating fleld is (Kul- 
srud and Pearce 1969) As ^ rL{E){B/5Bf, where rL{E) = E/ZeB = UlE2o.^/ZB^ Mpc, 
and 5B'^ = {kSBl)krr^=i. Here E20.5 = E/W^-^ eV and = E/IQ-^ Gauss. If \b{E) < koss 
(see expression 2), particles come only from a distance Dscatter{E) ~ V^shoss < hoss- This 
diffusive reduction in the observable volume can occur at the highest energies, including 
those above the saturation energy of the pion losses Eq — 2.S x 10^° eV, only if the general 
intergalactic magnetic field is strong, B > 2 x 10~^ Gauss, and the magnetic fluctuations 
are intense, SB/B ~ 1, conditions perhaps unlikely for a general intergalactic field. 

But at lower energies, the volume contributing to the observed particle fiux shrinks when 
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scattering is significant, compared to that observed if scattering is negligible. This shrinkage 
depends on particle energy, which means that diffusive transport can have a significant 
effect on the observed particle spectrum between the ankle at Ea = Q EeV and Eb = 
1.1 X 10"^^ {ZBgY^^ {6 B / B)"^/^ eV. The shape of the spectrum can be readily estimated in 
such a diffusive transport regime. The distance D{E) from sources to the observer is -De// = 
{XB{E)cTiossy^'^ = cTeff < cTioss- This contraction of the observable volume leads to the 
spectrum in (32) becoming J{E) fa [c/ ATx)Tf,f fUgUmdNi/ d{miC^'^) cx instead of oc 

E"^ ^ E^J oc in the energy range where pion losses dominate, E > A x 10^^ eV, 

instead of oc contrary to the observations. Solving the diffusion equation with systematic 
energy losses in the short mean free path regime confirms this estimate. 

Therefore, if there is a general intergalactic magnetic field, this and all other metagalac- 
tic source models which explain the absence of the traditional GZK cutoff by the existence 
of relatively nearby hard spectrum particle sources can work only if Eb < Ea, which im- 
plies ZBg{SB/Bf < 1.5 X 10"^ [equivalently, SB/B < 0.01/(^^9)^/2] _ ^^^^ler the general 
intergalactic field is weak, or the amplitude of the magnetic turbulence is small, or both. 

Compound diffusion, in which particles scatter strongly in localized regions, with scatter 
free transport between the scattering sites, yields a similar contradiction with the data. If 

the average distance between scattering centers is an energy independent distance / and the 
particles lose negligible energy while trapped in the scattering sites, then -De// = (IcTioss)^^'^ 
and and one finds E^JccE between Ea and the highest energy that can be trapped in a 
scattering site, also contrary to observations of the well-determined spectrum at 10 to 50 
EeV. 

The most likely candidates for such localized scattering sites are galaxy clusters, with 
magnetic fields perhaps a few xlO^'^ Gauss if smoothly distributed in the clusters (sec Eilck 
and Owen 2002 and references therein). The particles which can be diffusively trapped in 
a (large) cluster have energy less than 10"^^ ZB^^SB /By. Since E^J is fiat above Ea, one 
can conclude that if UHE cosmic rays arise in normal galaxies with active star formation, 
as is the case for the model proposed here, most of the particles get from the sources to 
us without becoming diffusively trapped in large galaxy clusters - which is not surprising, 
since galaxies with active star formation are not in rich clusters themselves, and rich clusters 
occupy only a small fraction of the metagalaxy. 

4. Acceleration and Escape of High Energy Particles from the Magnetar 

Wind and the Magnetar Wind Nebula 

If magnetars were born naked, the theory would be complete, in terms of its energetics. 
The model's basic results appear in Figure 1 and Table 1. For fiducial values of the input 
parameters and weak or moderate gravitational wave losses, the theory does too well - the 
source strength Kq is too large, leading to the inference that only 5-10% of the magnetars 
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are born with voltages sufficient to contribute to tfie UHECR. Tlie aim of tliis section is to 
estimate tfie losses imposed on the winds as they escape their natal supernovae, to address the 
acceleration mechanism, and to argue that only half the magnetars have magnetic geometry 
appropriate to their being UHECR sources. The losses associated with breakout from the 
supernova (estimated below to be ~ 30% of the initial rotation energy) affect the energy of 
the upper cutoff of the spectrum; the geometric factors affect the overall normalization. I 
also argue that in the equatorial outflow geometry appropriate to the model, adiabatic and 
radiation losses are negligible. 

4.1. Escape of the Wind from the Supernova 

Magnetars arc neutron stars. Neutron stars form in core collapse supernovae (Filippcnko 
2001 and references therein). The injection spectrum (23) can form only if the relativistic 
wind survives with negligible energy loss to large radii, and only if the high energy particles 
accelerated in the wind can escape into interstellar and intergalactic space with negligible 
adiabatic and radiation losses. In this and the following two sections, I outline the physics 
of the wind and particle escape which justifies adoption of the bare magnetar spectrum (32) 
as a reasonable "bottoms-up" model for the origin of UHE cosmic rays. 

Particles with energy E ~ E^j^^^^ an energy well above existing observations, must 
come from the outer limits of the relativstic wind, which at time t lie between r — ct and 
r — c{t — tspin,i), with tspin,i = niin(r£;M, tqw) - the outermost 5 light minutes or less. If the 
energy in this region instead is dissipated through doing work on the plasma and radiation 
in the supernova and the presupernova environment, the spectrum would be truncated at 
the highest energies, and possibly suppressed at lower energies if the environment dissipates 
wind energy at later times. The observed UHECR, with energies E < 10^°'^-E'2o.5 cV, require 
escape of the nonradiative wind to large radii only for Q < 3 x 10^ (£^20.5/^^1/^33)^^^ s~^; 
from (30), such rotation rates occur 45 minutes and more after the magnetar's formation. 

The essential issue arises in the initial spindown time being short or comparable to 
the time for the explosive ejection of the supernova envelope^. Magnetars, if born rapidly 
rotating, necessarily get rid of much of their rotational energy before the extended envelope 
of an isolated Type II SNe can be fully ejected, which opens the possibihty of the energy of 
the relativistic wind being dissipated in the dense matter around the newly formed magnetar, 
rather than being made available to relativistic expansion and particle acceleration at large 
radii. 



^Many of the topics discussed in this section have also been considered by Inouc et al. (2002) using a 
rapidly rotating pulsar with a 10^'^ Gauss magnetic field, in the context of the supranova model for gamma 
ray bursts, assumptions which yield rotational generation of electromagnetic energy on a time scale long 
compared to the supernova explosion itself. 
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If core collapse SNe form in binaries, or shed their extended envelopes before the core 
collapses through large mass loss, the presupernova star may be the compact {R^ ~ 10^'^i?5.3 
km) stripped helium core of the pre-supernova massive star, with envelope mass after core 
collapse M ~ 4M4 solar masses (Wheeler et al. 2000, Mazzali et al. 2002). Compact 
progenitors for core collapse SNe occur in Type lb and Ic supernovae, and it is interesting 
to note that the poorly known Ib/c supernova rates (Caparello et al. 1999) are comparable 
to the even more poorly known magnetar birth rate (Gaensler et al. 2001). Then the 
dynamical time of the envelope surrounding the newly formed neutron star is fairly short, 
td ~ 2.3RI^^ /mI^"^ minutes. The gravitational binding energy of the envelope is small, 
Wg ~ 2 X 10^°M|/i?5.3 ergs, a figure which includes the contribution from the neutron 
star's gravity. If there were no macroscopic electromagnetic energy input from the newly 
formed magnetar, the supernova shock, generated either from core bounce or neutrino driven 
convection and expanding through the envelope at speed Vg, would begin ejection of the 
envelope in tsN = R*/vs ~ I8-R5. 3/^530 seconds ^ td, where ^^30 = fs/30,000 km/s. The 
envelope thus forms a moving cavity containing the magnetar, expanding with energy ~ 
4 X 10^27^^4^^s3o ergs (Wheeler et al. 2000). 

The newly formed magnetar initially emits electromagnetic energy at a prodigious rate 
- Eem ~ 2 X lO^'^ ii'^^fl'l ergs/s, declining to ~ 10"^^ ergs/s after a few days. Thus, within 
the time used by the supernova shock to set the envelope in motion, the spinning magnet 
in the middle fills the cavity with "lightweight" relativistic energy, mostly in the form of 
electromagnetic fields, since at this early phase radiation losses bleed off the energy of any 
accelerated particles^. FiUing the cavity with buoyant "hghtweight" energy leads to prompt 
"shredding" of the envelope through Rayleigh- Taylor instabilities. 

Such shredding happens quickly. Consider the "radiation pressure" exerted on the en- 
velope by the electromagnetic fields spun off from the newly formed magnetar. This pressure 
creates an effective gravity at the expanding envelope's inner surface g = UEM-^enveiape/ ^envelope- 
Since this acceleration of the envelope's inner boundary is Rayleigh- Taylor unstable, buoy- 
ant "bubbles" of light weight energy rise through the shell at the generalized Alfven speed 
va = {Uem / PenveiopeY^'^ ~ [^EEAiit) / M^^"^ (adapting the results of Davies and Taylor 1950 
and Layzer 1954) and burst out, shredding the envelope in the time tghred- The total elec- 
tromagnetic energy contained after spin down to angular velocity is 

^Eem = Imlln l^^j ^ Imllnm ^ 1 X W^'^^^ ergs. (34) 

For n ~ ng{e2 = 1),E lO^o-^ eV, In A ^ 1.7, while at late times {Q < l^^, ^ < lO^^-^ eV) 
InA !=s lnA(n = 0) ~ 2.4. Expression (30) gives an adequate approximation for the time 



^Thc photon energy created contributes to the lightweight energy inside the cavity, until the Rayleigh- 
Taylor instabilities which shred the envelope allow all the photons, electromagnetic fields and relativistic 
particles to escape into the surrounding circumstellar/interstellar (ISM/CSM) medium. 
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to spin down to angular velocity Q; the general expression is in (26). Interstellar medium 
constraints (see §5) suggest gravitational radiation losses are substantial, thus requiring 
Qg < fli and the use of expression (34) in evaluating the electromagnetic input to the 
supernova. 

Rayleigh-Taylor instabihty, an adiabatic process, typically focuses the heavy fluid (the 
material in the envelope) into elongated "spikes" presenting a cross sectional size ~ R/k 
with respect to the outflow, covering a fraction ~ 1 /k^ of the initially unstable surface. The 
dominant scale of the magnetic Rayleigh-Taylor instability (Kruskal and Schwarschild 1954, 
Hester et al. 1996) as it becomes nonlinear is the largest unstable mode that will fit in 
the thick envelope, Aj ~ R/2. Therefore, the number of envelope fragments is N fragment ^ 
^envelope/ T^^l ~ 16. Electromagnetic compression rapidly collapses the fragments into a 
small percentage of the volume - simulations of magnetic Rayleigh-Taylor instabilities {e.g., 
Jun et al. 1995) suggest the fragments would have typical dimensions ~ 0.5Ai (k ~ 4) . 

Envelope shredding happens when the bubbles break through the expanding stellar en- 
velope. The bubbles rise faster in the expanding envelope as the electromagnetic energy accu- 
mulates. Most of the envelope's material is in a shell of thickness AR not greater than the pre- 
supernova star's radius i?*. Therefore tshred ~ AR/vA{tshred) ~ R*/vA{t shred)- Assume Tg ^ 
tshred > ^GW 7 assumptiou found to be self-consistent a posteriori. Then In A ^ 1 in (34), 
AEEMitshred) ~ (l/2)/02i^^ _ 4 5 ^ lO^^ (//i3//45e^) lu A, va « 7500(//33/e2) A//45M4 
km/sec and 

— ) seconds; (35) 
InA J 

this estimate assumes the envelope thickness is comparble to its radius. 

The energy AEghred expended in shredding the shell into fragments which cover a frac- 
tion of the area around the newly forming relativistic bubble is not large. The mechanical 
work that has to be done to shred the shell is 

AE shred ~ Mgnvelope9'^A{tshred)t shred ~ AEem 

Vs 

= AEem {AEem/MvIY'^ \t=t,^^^, = (1 - W shred) AEEM\t=oo. (36) 

I used g — UemA/M ^ AEEM/MR{t) ^ AEEM/Mvst to obtain this expression. 

Shredding of the envelope leaves most of the electromagnetic energy free to escape into 
the surrounding interstellar medium. From expression (36) the fraction 

w - 1 - ^^'h^^d _ 1 m 2.4 / /^L \ ^'^ .oyx 
AEEM\t=o. '•'lnA(Q = 0) U45M46i^ioy ^ ^ 



of the electromagnetic energy can escape to form a relativstic wind and Magnetar Wind Neb- 
ula, assuming mixing of the nonrelativistic material from the shredded supernova envelope 
does not create too much mass loading. 
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The relativistic fields and particles emitted by the magnetar behave as a relativstic 
"fluid" only so long as the non-relativistic material in the envelope remains unmixed with 
the relativistic gas. The Rayleigh- Taylor "fingers" are themselves subject to disruption 
from subsequent Kelvin-Helmholtz instabiltities, whose final consequence is to mix the non- 
relativistic material of the presupernova star's envelope with the relativistic outfiow im- 
mediately in contact with the nonrelativistic material, creating local "mass loading" . The 
Kelvin-Helmholtz mixing speed when one of the fluids is relativistic is (Arons and Lea 

1980) Vmi,^ ~ O.ljKH/k ^ Q.liUEM/Pfingerf'^ = ■IVa/k^''', whcrC Va = {AEem/MY/'' 

as in the discussion of interchange motions. For the environment considered here, the mix- 
ing times are on the order of an hour, during which time the envelope fragments move 

to several hundred times the size of the presupernova star, while the relativistic gas not 
in contact with the fragments escapes at the speed of light. The nonrelativistic matter 
mixes into only a small fraction, ~ ~ 10-20%, of the relativistic outfiow stream- 

lines, thus leaving almost all of the relativistic energy that survives envelope shredding 
and fed in by the magnetar at later times to be free of mass loading and able to acceler- 
ate high energy particles at radii large compared to the size of the ejecta. The fraction 
of the energy surviving mass loading thus is given by the covering factor of the shredded 

ejecta, Wioad ~ 1 - Nfragmenti^ fragment /^envelope) = 1 " {N fragment /^^k"^) ~ 0.75. ThuS, 

Wbiowout = WshredWioad ~ 0.6 - 0.7; 60-70% of the lightweight energy survives to break 
out of the confining supernova to blow a cavity in the surrounding interstellar/circumstellar 
medium (ISM/CSM), before the envelope ejecta get a chance to drive the standard super- 
vova blast wave. Instead, the relativistic bubble (large scale, low frequency EM fields and 
particles) drives a relativistic blast wave into the ISM/CSM. 

The losses involved in the relativstic wind pushing out of the supernova come from the 
initial rotational energy loss put into the wind. Therefore, the highest voltage parts of the 
wind would be lost, which would reduce the maximum particle energy that can be achieved 
by a factor W^H^^, ~ 0.8. 



4.2. Relativistic Blast Wave and Magnetar Wind Nebula in the Interstellar 

Medium 

The particles accelerated in the wind still must escape the Magnetar Wind Nebula 
without appreciable loss, even after the wind and the MWN escape the supernova, if the 
bare magnetar spectrum (23) is to represent the rate ultra-high energy ions are injected 
into the metagalaxy. For times greater than tshred, the relativistic energy escaping from 
the supernova drives a shock wave directly into the surrounding interstellar or circumstellar 
medium, superseding the usual shock driven by the supernova ejecta - the ejecta are left 
behind. For simplicity, I consider only the case of an explosion into a uniform external 
medium of mass density pi, the results for an explosion into a cavity excavated by mass loss 
from the supernova's progenitor, with mass density oc r~", are similar to the uniform case. 
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Since AEEM/{ct)^ ^ pic^ for many months after the explosion, for all reasonable 
densities of the interstellar medium, the relativistic energy drives a relativistic blast wave, 
whose Lorentz factor and radius are 



Rh = 



17 AEem 
87rpic2(ci)3 

t- I dt' 



1/2 

(38) 



2r^(f) 



ct 



(39) 



The dependence of F;, on AEem and pi comes from simple consideration of momentum 
conservation; the factor 17 /Sn comes from the similarity solution (Eltgroth 1972, Blandford 
and McKee 1976). The shocked interstellar medium occupies a thin shell, whose thickness is 
^Rism — Rb/Q^bit)- Within this thin shell, the interstellar magnetic field is greatly amplified, 
B2ism — 3i?ir(t). Flux conservation since the start of the blast yields the coeffcient 3 in this 
expression; the instantaneous jump conditions replace 3 with 2, which applies to the B field 
just behind the blast wave. The relativstic expansion continues until the time Trei such that 
^b{Trei) ~ 1, which yields 



17 AEemV^"^ .^^fAE; 



rii is the baryon number density upstream of the blast wave in units of 1 cm~^. For t > T^eu 
the blast enters the nonrelativistically expanding, energy conserving Sedov phase. 

The relativistic fields and particles from the magnetar fill the volume out to a contact 
surface located at Rc = Rb — ARism, which is close to Rb when t < Tj-ei', it suffices to 
take Rc as being equal to Rb- This volume forms the Magnetar Wind Nebula (MWN) 
in its relativstic expansion phase. The freely expanding magnetar's wind terminates in a 
relativistic reverse shock wave located at radius Rg where the momentum flux in the wind 
matches the pressure in the magnetar wind nebula (almost all of which is in the initial energy 
pulse WbiowoutAEEM) ■ I consider only times greater than tshred > ^Gw, since at earlier times 
the relativistic blast wave has not yet formed. 

The pressure balance condition is 

EEMit-Rs/c) 3 \mlWblo.nout\^^m 

One readily finds, upon using (30) and assuming Rg is not almost equal to ct, that the reverse 
shock forms where 

RA' 2 /2r,^^^-^ 



1/3 



(42) 

Rb J SWblowout^''^ \ t 

with p = 1/4 for tqw -C ^ < 2Tg, and p = 1/2 for t > 2Tg. At early times, when gravity 
waves control the spindown, the reverse shock expands with the MWN (4p —1 = 0, Rs — 
{2\^Wbiowout In A./3)ct 0.63ci), but for times greater than 2Tg, p — 1/2 and the blast wave 
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leaves the reverse shock behind^, with 

Because the wind, at r < Rs, and the MWN, between the Rg and Rc, both expand rela- 
tivistically, the electromagnetic fields everywhere (0 < r < Rc) have amplitude proportional 
to 1/r, 

where a is the ratio of the electromagnetic energy flux to the plasma kinetic energy flux in 
the wind interior to Rg. f{r,t) describes short wavelength [A ~ 27ic/Q(t — r /c)] structure in 
the wind, as is the case if the wind has wave structure similar to vacuum strong waves or, 
if an MHD wind, is "striped" (Michel 1971, Coroniti 1990, Michel 1994, Bogovalov 1999), 
with the B field direction reversing with wavelength A while the overall amplitude declines 
inversely with r. 

Pulsar Wind Nebulae (PWN) observed in our own galaxy expand nonrelativistically, 
due to the fact that the initially injected relativistic energy AEem is small compared to 
the kinetic energy of the conventional supernova. The resulting "bags" of electromagnetic 
energy and particles injected by the pulsar are confined by the inertia of the supernova 
ejecta until times large compared to their ages when we see them, as is the case in the Crab 
Nebula (Hester et al. 1996), or by the surrounding CSM or ISM, as might be the case for 
the nebula around PSR 1509-58 {e.g.^ Gaenslcr et al. 2002). In contrast, the hypothesized 
MWN described here expands relativistically at all radii Rg < r < Rc = Rb ~ ARism, 
for t < Trei- Therefore, the post-shock velocity at Rg must be relativistic, and a ^ 1 - 
the magnetar's wind and the MWN must be everywhere electromagnetically dominated and 
"clean" - throughout the relativstic expansion phase, the field structure will not be disturbed 
by interaction with clumps of nonrelativstic matter, since all such material lies in the shreds 
of supernova material through which the wind escapes, and in the shocked shell of ISM/CSM 
right behind the blast wave. 

If magnetars behave like pulsars and have wind outfiows with a dense plasma {N± ^ 
Ngj) in addition to the Goldreich-Juhan fiuxes associated with the electric currents (the 
high energy ions form one of these electric current fiows), the wind and MWN would have 
a structure perhaps described by MHD. N± for a magnetar is not well understood, either 
theoretically or observationally. Observed magnetars are slow rotators with voltages below 



^If gravity wave losses arc negligible, pure electromagnetic spindown yields {Rg/Rbf' = 
{2liWuowout){2TEMlt), Rs = V^/^Wbiowoutc{2TEMt)^/^ = 4 X lO^^hs / ^ih^lV^Xilrs Cm when 
t > 2tem > tshred- Since tshred < 2tem, at very early times Rs « {2/3Wbiowout)^^'^ct = 3 x 
10i3(o.7/Wwo^o„t)V2(t/iOmin) cm, t < 2tem. 
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the level at which pulsars stimulate observable PWN, so prominent nonrelativistic MWN 
around observed galactic magnetars would not be expected. Theoretically, pair production 
in the ultrastrong magnetic fields of magnetars has been controversial, with the inhibitions of 
pair creation by photon splitting being uncertain (see Baring and Harding 2001). However, 
since a must be large, the MWN would not be an efficient radiator, if all the acceleration of 
pairs is associated with the reverse shock at until t > T^ei ~ 1 year. 



4.3. Acceleration of UHE Ions in the Magnetar Wind 

The metagalactic magnetar model for UHE cosmic rays has been developed using the 
phenomenology and theory of relativistic winds in galactic Pulsar Wind Nebulae (PWN) 
as the basic physics input. So far 1 have not described a specific acceleration mechanism, 
through which the ions in the equatorial Goldreich-Julian current actually tap the voltage 
available in the wind; the same lacuna exists in the theory of ions in PWN. A detailed 
investigation of this subject is outside the scope of this paper. However, a few remarks are 
relevant, addressed to the most likely mechanism in this context, the "surf-riding" of charged 
particles in the relativistically outflowing electromagnetic flelds. 

As is clear from expression (7), the ions would gain energy comparable to the total 
available voltage if they can cross magnetic fleld lines for a distance comparable to one 
decade in radius. Such fleld line crossing can occur in either radius or in latitude, or both. 
To see how the simplest form of surf-riding works, consider particle motion in the asymptotic 
wind of the aligned rotator (Buckley 1977, Contopoulos and Kazanas 2002). 

Substantial theoretical evidence has accumulated that when the electromagnetic flelds 
dominate the outflow's energy density, the structure of the dipolar magnetosphere and its 
outflow can be approximated as force-free [pgE+ (l/c)J x B = 0, with pq the charge density 
and J the current density], and that under the magnctohydrodynamic conditions appropriate 
to PWN {n ^ Pq/<li where n is the total quasi- neutral plasma density), the structure of the 
aligned rotator's wind at radii larger than a few closely approximates that of the split 
monopole (Michel 1974 , Contopoulos, Kazanas and Fendt 1999), except perhaps in narrow 
cones around the rotation axis. These flelds have the form (Michel 1973, Bogovalov 1997) 







B4 


MsmO 






Eg 


— B^ 






Be 


— Ej- = 


Ea, 


= 0. 



The magnitude of the monopole moment is M = kfi/R^, where p is the dipole moment and 
A; is a constant on the order of unity. For the self-consistent numerical solution obtained 
by Contopoulos et al. (1999), k = 1.36. The coordinate system is spherical, with 6 = 
being the colatitude of the angular momentum vector, and with the azimuth being 
measured with respect to an arbitrary x axis. The plus and minus signs apply to the opposite 
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( northern/southern or southern/northern, with respect to the rotation axis) hemispheres. 
The fields reverse direction abruptly across the rotational equator {9 — 7r/2), which implies 
the existence of an equatorial current sheet where the return current flows. In the body of 
the wind, the volume current is radial. 



The electromagnetic energy propagates with the fleld fine velocity 



E X B I ^ x" X 

Ve = C = C ( T— + 
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(46) 
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(47) 



Charged particles flow out with with velocity v — VE + v\\b, where b is the unit vector along 
the magnetic field. When particles enter the wind with initial energy mc^^L (corresponding 
to initial flow 4- velocity c(3l^L: 7l = (1 — /?1)~^^^), their parallel speed is, when 7l ^ 1 
(Arons, in preparation) 



7l 



(48) 



Therefore, the wind's Lorentz factor is 

7. = (l-/5|-/3|r^/' = (7i + 7l)^/^ 



7l 



— 1 sin^ 9 

R-L 



1/2 



(49) 



(Contopoulos and Kazanas 2002). When rsin^ 'JlRl: 7 — rsm9/RL (Buckley 1977); 
the wind is a linear accelerator, with each particle's energy increasing by one factor of 10 for 
each decade in radius. At large radii, the particles move with the fleld lines, with negligible 
motion parallel to B, that is, they "surf-ride" on the electromagnetic fleld. 

This result applies to quasi-neutral winds, whose total density is large compared to the 
Goldreich- Julian density; it does not apply to the charge separated force free split monopole 
(Arons, in preparation). Pair creation within the magnetosphere supplies densities large 
compared to ugj, with injection energies 7^ ~ 10 — 10^; the precise values of n and '-/l 
depend on the specifics of the magnetospheric parameters (surface magnetic field strength, 
voltage and whether curvature or inverse Compton gamma rays radiated by electrons or 
positrons accelerated within the magnetosphere initiate the pair cascades.) 

Such surf-riding acceleration may seem surprising, since v ■ E = 0. The physics of the 
acceleration appears when one considers the (small) effect of inertia on the velocity transverse 
to B, by solving the equation of motion to first order in the rest mass m. The solution for 
the velocity perpendicular to B at large radius, where the velocity becomes equal to the 
radial component of the E X B drift, is 



V 



VEr + 



mcB X D{'jv)/Dt 



VEr + 



mcB X D{jEVEr)/Dt 
~~Q 52 



= VEr + Vpol- 



(50) 
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Here D / Dt — d/dt+v ■ V ^ d/dt+VEr{d/dr). The polarization drift velocity Vpoi is parallel 
to the electric field: 

v^oi = — ^-B X — — = =F ^, (51) 

qB^ Dt ^ q M ' ^ ' 

where 6 is the unit vector in the 9 direction, and I have assumed rsin^ ^ R^. Vpoi always 
has a direction such that qVp^i ■ E is positive - particles gain energy from the electric field in 
an accelerating fiow. Then energy conservation takes the form 

2-^7 ^ mc^ sin 6* . ^. 

mc — = qvpoi ■ Ee = H . (52) 

Dt Kl 

Integrating with D/Dt = c{d/dr) yields expression (49) in the limit rsin^ ^ '^lRl- This 
correspondence can be extended to all radii, at the expense of much more complex algebra. 

Thus relativistic force free winds can be linear accelerators, with the particles "surf- 
riding" on the electromagnetic fields. For the aligned rotator, these fields form a (spherical) 
DC transmission line, with the fluid accelerating following the increase of the energy trans- 
mission speed toward c as the magnetic field approaches a purely transverse structure with 
increasing radius - the wind 4-velocity 7^ is not constant but is asymptotially oc r. For 
c/VL = 3O/Q4 km and 7^ < 10'^, such linear acceleration sets in at radii not larger than 
(30, 000/1^4) (7l/10^) km. 

The result (49) applies to force free conditions, which require a = B"^ /Anp^fyjC^ ^ 1, 
where p is the rest mass density measured in the frame where the neutron star's center 
of mass is at rest. Assume plasma at nonzero latitude |A| = \9 — 7i/2\ > crosses the 
light cylinder with initial energy /particle uic^^l — El and rest mass density Pl — {iTii + 
2K±m±)nGj{RL) = meff^B{Ri,)/2Trcq ~ meff^n/2T:cqR\, where I have assumed the heavy 
ions have rest mass rrii and density rij = noj and are mixed with a pair plasma of density 
n^-\-n- — 2K±naj ^ n^. Then 

^ ^ 47rme//(Q-Bi/27rc5)7Lc2 2meff(rfL^ 2me{fC^^L 

rueff = mi-\-2K±m±. (54) 



In the absence of magnetic disspation and for rsin^ ^ ^lc/VI, a — (TL'^L{RL/rsm.6), 
with the force free region ending aXr = Ri — Rl'^l'^l/ sin^, where a drops to unity. Beyond 
Ri, the plasma probably coasts, with 7tt, = 71 ~ ^Lf^L-i in regions where sin^ ~ 1. If the ions 
are absent, me// — >■ 2m±K±. 

The polarization drift velocity (51) causes the ions (and other particles) to drift across B 
(that drift is the origin of the acceleration), moving off their radial orbits toward the rotation 
pole a distance A/(r) = {mrestc^ / (l^mag){''^^ / Rl)- Since the acceleration ends at radii no 
larger than r = Ri, the maximum fractional energy gain for an ion from polarization drift 
alone is E/q^^^aag = mi/rrieff- If mg// ~ 2 — lOmj and 7^ ~ a few (as is the case when pair 
creation cascades are strong, as in the young pulsars) in the return current region, 10% - 50% 
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of the magnetsopheric potential on open field lines could go into accelerating the ions due 
to polarization drift along E. This fraction is lower if the pressure of the MWN terminates 
the wind interior to Ri. 

The accelerator closely approximates a purely radial motion - curvature (a.k.a. syn- 
chrotron) radiation losses are neghgible, so long as the particles surf- ride {v ^ vet)- 

The energy that can be achieved within r — Ri thus depends on the plasma composition. 
For an aligned rotator with magnetic moment parallel to the angular momentum, the ions 
in the wind would flow in the equatorial current sheet, while at higher latitudes, the plasma 
would be composed solely of pairs. The theory of pair creation in ultra-magnetized neutron 
stars is not as well understood as in standard pulsars ~ 10'^'^ cgs), since a variety of 
strong field QED effects can intervene to inhibit pair creation. A traditional estimate of k,^ 
which ignores these effects {e.g., Ruderman and Sutherland 1975, Arons and Scharlemann 
1979) and that takes account of the limitations set by curvature gamma ray emission on the 
energies of electrons in the primary polar particle beam yields k± ~ lO'^di^^Rl^y^'^. Taking 
account of photon splitting, and pair production directly into the lowest Landau levels, 
effects that reduce the contributions of the synchrotron cascade to pair creation (Baring and 
Harding 2001), suggests k± ~ 10 — 100. This polar flow would reach the light cyhnder with 
Lorentz factor ~ 10 — 100, yielding Ei = m±c^7i = m±c^'jL(^L = 3 x 10'^^ fi^sVlKlO / k± for the 
asymptotic 4 velocity of an ideal MHD, pair dominated wind outside the equatorial current 
sheet, assuming no magnetic dissipation interior to Ri. 

In the aligned rotator models contructed to date, the dynamics of charges in the current 
sheet, where the ions flow as the return current, are not specified. Simple models of the 
current sheet's structure (Arons, in preparation) suggest it to be a force-free rotational 
discontinuity separating the hemispheres. This structure causes the ions to be frozen to the 
fields which partcipate in the general linear acceleration, and therefore also obey the linear 
acceleration law (49). If the geometry has the dipole moment antiparallel to the rotation axis, 
then the ions are themselves the charge carriers in the |A| > hemispheres and participate 
in the generaL acceleration, while the current sheet carries a negatively charged outflowing 
return current. 

In general, magnetized rotators are not expected to have their dipole moments exactly 
aligned or counter- aligned to the angular velocity. While a solution of the force-free oblique 
rotator with dipolar field structure inside the light cylinder has not yet appeared, the ex- 
perience with the aligned rotator suggests the wind emerging from an oblique rotator will 
have the structure of the oblique split monopole at radii larger than a few times Rl- An 
exact solution for the oblique split monopole has been obtained by Bogovalov (1999), whose 
results for the fields are shown in Figure 2. 

The magnetic and electric field strengths in the oblique split monopole are the same as 
in the aligned split monopole; the only difference is that in the equatorial region (latitudes 
— i < A < i, where i is the angle between the magnetic and rotation axes), the directions of 
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E and B reverse every half- wavelength, where the wavelength is 2'kRl, forming a "striped" 
wind (Michel 1971, Coroniti 1990). The current sheet is now wrinkled and is frozen into the 
flow, advecting outwards with the same 4-velocity as the general wind. That 4-velocity is 
exactly the same as in the aligned split monopole, given by (49). Thus, under the plausible 
assumption that the obhque dipole's wind takes the asymptotic form of the oblique split 
monopole, the maximum ion energy/particle achieved in all magnetic geometries of rapidly 
rotating magnetars through the action of polarization drift of the ions along the MHD electric 
field is 

Emax = miC^-fi = nXiC^-fLfyL = q^mag —, (55) 

^eff 

where mgff = rrii + 2n±m± and n± is the pair multiplicity in the current sheet. 




Fig. 2. — Frozen- in current sheet structure of the oblique split monopole in the inner wind, 
from Bogovalov (1999). Left panel: Meridional cross-section of the poloidal field structure 
at large r, showing the crinkled current sheet. Right panel: Intersection of the curent sheet 
with the equatorial plane. The toroidal magnetic field forms stripes with opposite directions 
between each current layer, as indicated by the arrows between the current sheets. 

If K-i- < 10'^ in the current sheet, m ~ mj and this elementary form of surf-riding 
yields aL = 0.5Ze<!>mag/Ampc''^L = 6 x 10^\Z/A)^i33^l{10/^L). Then R, = RlIlcJl = 
3.5{Z/A)fissnl pc, and Emax = 3.3 x lO^^Zfi^snj eV. 

Linear acceleration due to polarization drift cannot not go to completion in the freely 
expanding (r < Rg) magnetar wind since Rg -C -Ri. The MWN exterior to the large a 
shock ending the wind also expands relativistically but with a decelerating field line velocity. 
Therefore, if the magnetar wind model is to explain UHECR, there must be a still stronger 
acceleration gradient in the wind; laminar surf-riding in the decelerating MWN saps particles 
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of their energy^°. 

The polarization drift based hnear accelerator is only a lower limit to the rate at which 
surf-riding in the magnetar's wind can can accelerate the charges. If the fields become 
dissipative within the wind, more rapid acceleration can occur. The low magnetization 
inferred for the equatorial sectors of pulsars' winds suggests that such dissipation must 
exist. 

Figure 2 shows that for an obique rotator, much of the Poynting flux is tied up in 
waves. In the ideal MHD outflow, these waves are frozen in current sheets. These structures 
become thinner as the wind expands, and are subject to inductive dissipation (Coroniti 1990, 
Michel 1994, Lyubarsky and Kirk 2001, Lyutikov 2002) at and beyond the disipation radius 
of the crinkled current sheet Rdiss ~ '^i^±Rl ^ Rs- They are also subject to instability with 
respect to emission of kinetic Alfven waves (Bellan 1999, 2001, and Arons, in preparation) 
which propagate with respect to the plasma in the wind frame. Current driven instability 
with respect to wave emission is a faster process than the inductive dissipation mechanisms 
discussed to date. 

Assume wave emission is the dominant process dissipating the wrinkles in the current 
sheet, and that such emission leads to the formation of large amplitude electromagnetic 
waves with dimensionless amphtude a — qSB/rriicQ » 1. These waves can accelerate ions 
faster than does polarization drift in Eg. Even when r < Rdiss, the Alfven waves can act 
as powerful ponderomotive accelertors (Chen et al. 2002), while for r > Rdiss these waves 
become superluminal "strong waves" (Melatos 1998), whose ponderomotive force also yields 
surf-riding acceleration. When the waves are kinetic Alfven modes, the ponderomotive force 
applies preferentially to the ions. 

An estimate of the ponderomotive work done suffices for the present discussion. 

This single particle force arises since the transverse motion of a charged particle in the 
wave electric field gives rise to v x B force in the wave magnetic field which points radially 
outwards, does not average to zero over a wave oscillation, and has magnitude 

fpond = -{Sv x5B) = ^ ^sb) = ^ (— X Sb) = mcn^£=r, (56) 

c mc \-fosc I mdl \-fosc / V 1 + a 

where a particle's transverse oscillation momentum in the strong wave is Sp^ ~ qSE/Q, the 
oscillation energy per particle is niiC^'yosc = miC^[l + ((5pj_/mc)^)]^/^ = rriiC^y/l + a^, r is 
the unit vector in the radial direction and ( ) indicates an average over wave phase. The 
work done on a particle as it moves a distance Ar through the wave is then, when a ^ 1, 

SB Ar Ar 
Wpond = f ponder = {mcQ a)Ar = qBr— ^ q^magV, V = — , (57) 



^°In pulsar wind nebulae, large-scale surf-riding in the wind is an adequate explanation of the infered wind 
four velocities (Contopoulos and Kazanas 2002, Spitkovsky and Arons 2003) 
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since 5E ^ 5B and 5B/B > 1 in the rotational equator - most of the energy flow is in the 
wave Poynting flux. The evidence from Pulsar Wind Nebulae is that 7/ ~ 0.1 and therefore 
Ar ~ O.lr. The fact that Ar <^ r and the particles gain only a relatively small part of 
the total voltage available most likely is due to the demise of the strong waves coherence, 
as a result of parametric instabilities in these electromagnetic structures - see, for example, 
Leboeuf et al. (1982). 

This argument suggests that such surf riding may be the source of the acceleration 
hypotheiszed in the model for UHECR studied here. A detailed demonstration that this 
is so is a task for another study. If so, the ion acceleration is prompt, beginning a few 
seconds and continuing for a few days after the escape of the relativistic wind from the natal 
supernova - that escape occurs ~ a minute after the neutron star's formation - see expression 
(35). 

The standard alternate for the acceleration mechanism of UHECR is some form of shock 
acceleration, usually invoked as diffusive Fermi acceleration in relativistic shocks behaving in 
a manner similar to nonrelativistic shocks {e.g., Waxman 1995, Bahcall and Waxman 2003); 
see Ostrowski and Bednarz (2002) and references therein for a discussion and critique of 
this mechanism's plausibility in relativistic shocks.) In the current context, the fact that the 
electromagnetic fields dominate the energetics of the outflow militates against the relevance of 
that mechanism. Alternate, "flare" mechanisms, based on direct dissipation of the currents, 
are even more poorly understood. In this regard, identifying bursts of UHE cosmic rays 
as coming from the same direction as bursts of quasi-coherent gravitational radiation, with 
correlated arrival times, would give substantial empirical input into constraining the specific 
physical processes behind the acceleration of the highest energy cosmic rays. 



4.4. Escape of UHE ions from the Magnetar Wind Nebula 

In accord with the discussion in §4.3, I assume the high energy ions achieve their energy 
in the relativistic wind through surf-riding acceleration, at a radius /2„ ~ Rdiss ~ i^±Rl = 
(10^ - W)Rl = (30,000 - 300,000)7^4 km. Beyond Ra, dissipation of the strong waves 
reduces the structure of the magnetic field to that of the ahgned rotator (Blandford 2002), 
with fields given by (45) but now with M = ji / Rl ^ M\\ = ji cos i/Rl. The acceleration 
injects the ions into the angular sector —i<\<i, with a 4-velocity large compared to the 
wind's 4-velocity (49). One can readily show, either from the equation of motion of a UHE 
ion in the crossed E and B fields of the wind, or from the motion of an ion in the pure B 
field in the E — Q frame of the wind, that the radius of curvature of an accelerated ion's 
orbit is 

Pc = 2ri7^ = (2r/7^/ cosi)r, (58) 

where tl = rriiC^'-f/qB = rjr is an ion's formal Larmor radius. When fa r/i?; ^ 1, as it is 
at radii comparable to and larger than Ra, Pc ^ r. 



-34- 



This large radius of curvature of ions' orbits, applicable wherever the flow is relatisti- 
cally expanding (both within the wind and in the MWN) means that both synchrotron and 
adiabatic losses of the ions are negligible. From the power radiated by a charged particle 
with orbital radius of curvature pci P = (2g^c/p^)7^, 7^,, ~ r / Rl, 7 = rjq^/miC^, one readily 
finds 

^ 6 A'^ f m„c^\^ m„c^ ^ ■) . f r \^ 3.6 A^sec^i a . . ^, 

where rji = rj/O.l and Tiq = r/10^° cm; one also obtains this same radiation time by applying 
ordinary synchrotron losses in the fluid frame of the wind where E = 0. Thus, for r > Ra > 
10^° cm, radiation losses are negligible {Trad ^ r/c). Also, since p ^ r, the high energy 
ions are not coupled to the field lines and do not suffer adiabatic expansion losses, so long 
as they escape the wind and the MWN while the MWN's expansion is relativistic. 

The particles leave the wind and move through the MWN with speed v{t) = c(l — 
1/27^) pa c. For the observed UHE cosmic rays, with Ea < E < 2x 10^° eV <C E^ax, the par- 
ticles at radii large compared to Rg were accelerated at radius Ra {Ra < Rs by assumption) 
and at the time r„ = Ra/c+ {Ea/E)ta, where i„ = {9/8){Zer]/ pEa) ~ 2Zr}ih^/ p^^Ea,is.s 
days is the time at which the neutron star's voltage was at the value corresponding to ac- 
celeration of ions at the ankle in the spectrum. Here i?a,i8.8 = -E'a/6 x 10^^ eV. For particles 
in the UHECR spectrum measured to date, fa ta{Ea/E), so long as Ra < Rs- At time 
t, particles of energy E have reached Rparticie — c{t — Ta) + Ra] for all particles of energy 
E>Ea,i = E/mc^:>Tb{t). 

The outer boundary of the MWN lies at Rb{t) given by (39). All particles with energy 
E > Ea eV catch up with the blast wave at the time Tcatchup when Rparticie first exceeds 
Rb{t). Equating the radii yields 

Tcatchup = ^Ta =8.3 '■ — months. (60) 

V TT piC^ J \ ni P33Ea,18.8 E J 

Here AEem,52 = AEem/IO^^ ergs. 

Comparison of (60) to (40) shows that all the UHE cosmic rays reach the blast wave 
before the MWN becomes nonrelativistic. The 4 velocity of the blast wave at Tc is Pb^i,{Tc) — 
4:{E52/niy/^{p33Ea,i8.8/ Zr]ih5)^/^{E / Ea)^^^ . The cosmic rays have no difficulty getting through 
the plowed up interstellar medium and its compressed magnetic field just behind the blast 

wave, once they arrive at the MWN's outer boundary - the cosmic rays' Larmor radii in 
the shocked interstellar magnetic field are r^s = 5A0{Ea,is.8/Bi,^){A/Th[Tc{Ea)]){E/Ea)^/^ 
parsecs, which is far larger than the thickness of the shocked shell of interstellar matter. 
Thus the UHE ions avoid substantial adiabatic energy loss in escaping the MWN. Here 
Bi,^, = Bism/'^ /xGauss. 

Other radiation losses in the MWN might occur if the high energy ions encounter a 
sufficiently dense external radiation field. The supernova's own optical and ultraviolet lumi- 
nosity is the most prominent candidate for such radiation. A reasonable upper limit for that 
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luminosity is ~ 10^^ ergs/sec, emitted over the first few days, i.e., over the same time as the 
UHECR acceleration occurs in the wind. At the radii where the acceleration is likely to occur 
(see §4.3), photopion losses in the stellar radiation field might limit the maximum ion energy 
that can be achieved to no more than 10^^ eV. Losses due to Compton scattering and pair 
creation are negligible, as is energy loss in the dense matter of the supernova, since the ejecta 
have only a small filling factor, and move out slowly behind the high energy particles. Losses 
in the radiation field of the MWN are unhkely to be of significance, since this highly mag- 
netized flow is unlikely to be radiatively dissipative until after the outer blast wave becomes 
nonrelativistic, which occurs after the UHE ions have left the system. Radiation losses in 
the interstellar environment surrounding the MWN (pion and pair production, synchrotron 
radiation, inverse Compton scattering) are all negligible, even if the MWN forms within a 
cluster of hundreds of O and B stars, and the magnetic field strength is that of molecu- 
lar clouds (milligauss). Interactions with the circumstellar matter that might surround the 
presupcrnova star also can be neglected. 

Radiationless escape of the UHE ions from the MWN works only when the ions form 
the equatorial return current. This geometry applies to rotators with dipole axes making an 
acute angle with respect to the rotation axis. Thus, only 50% of the magnetars can supply 
UHE cosmic rays, contributing the factor Wgeom = 0.5 to the flux of UHECR (expression 
32). 

If the rclativistic wind has a dense pair plasma, as well as the ion return current in 
the equator and the polar electron current, Compton scattering from the pairs might create 
a high energy photon flash, which may be an interesting candidate for a prompt Gamma 
Ray Burst. This is a topic for a separate investigation. For the present purposes, it suffices 
to observe that the supernova photon luminosity is small compared to the magnetar's spin 
down luminosity, therefore Compton drag can have little effect in dissipating the electro- 
magnetically dominated outflow. 

5. Interstelletr Supershells and Supernova Remnants 

Formation of interstellar HI shells and supershells provide an interesting constraint 
pointing toward only a small fraction of the newly born magnetars having voltages large 
enough to accelerate UHECR, and perhaps to large gravitational wave emission by such 
objects. 

If the magnetars lost energy purely electromagnetically, they would deliver energy on 
the order of {l/2)IQ,^Wi,iou,out — 3.5 x 10^^ I 45Q,l{Wj,iou,out/^ -7) ergs with each event. These 
explosions would generate supershells, large expanding structures in the interstellar medium 
with energy in excess of 10^^ ergs (Heiles 1979, 1984) and dimensions in excess of 100 pc. 
From Chevalier's (1974) expression (26), with a random velocity of clouds in the ISM of 
~ 20 km/s, electromagnetic energy deposition with 1^4 ~ 1 would create shells with radii 
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~ 230 pc. Further investigation over the last 20 years has led to the conclusion that all 
such events in our galaxy occur at a rate I'supersheii ~ 0.8 x 10~^ yr~^ (Ehlerova and Palous 
1996), about a factor of 12 slower than the current estimate of the magnetar birth rate; 
other relations between explosion energy and remnant size {e.g. McKee and Cowie 1975) 
lead to similar conclusions. Shells as large as 200 pc form only a small, but ill determined, 
fraction of the total supershell population. Therefore, if the theory outlined here is to be a 
viable explanation of UHECR, either the magnetar rate is overestimated, or the energy per 
magnetar delivered to the interstellar medium is smaller (or both). 

The approximate correspondence shown in Figure 1 between the theory outlined in this 
paper and the observed UHECR spectrum requires = Wgeom^mAihb/ fJ^33)ng2THu ~ 0.03 
if the spindown is purely electromagnetic. This result comes from requiring the model to 
explain the well determined spectrum between 10^^ and 6x10^^ eV. If the particle spectrum 
measured at energies below the ankle extrapolates into the UHE regime, Kq drops to 0.02 for 
the pure electromagnetic case. If one restricts oneself to magnetars with parameters defined 
by current observations, the main freedom is in the fact that we don't know whether all 
such objects are born with the same initial rotation rate. Magnetars can succeed as a model 
for UHECR if their voltages are sufficiently large, requiring initial periods Pj < 4(771/^33)^/^ 
milliseconds for at least some of the magnetar population. 

With Wgeom = 0.5, and rig/0.02Mpc^'^ ~ 1 , the inferred value of Kq tells us that 
the birth rate of such rapidly rotating objects, which by hypothesis are subject only to 
electromagnetic torques, is z/m"**'' ~ O.lz/m = lO'^Um^ yr^^ (see §3.3). Then most of the 
supershells in our galaxy would be the result of these "electromagnetic bombs" in the inter- 
stellar medium, and most would be of the largest variety. Loeb and Perna (1998) made the 
related suggestion that Gamma Ray Bursts drive supershells. 

If one insists that most supershells do not come from exotic sources (probably true), and 
that their numbers not be dominated by the largest variety (certainly true), then the energy 
delivered to the interstellar medium in each event is less than the (1/2)1^2^ = 5 x 10^^/45^1 
ergs. For pure electromagnetic spindown, initial rotation rates less than 10^ sec~^ are the 
only means of achieving that reduction. Creating shells smaller than 100 pc, which would 
cause them to be missed in Heiles's (1984) census, requires the energy deposited per event 
to be less than 4 x 10^^ ergs, in turn requiring Qi < 3 x lO^sec"^. Then the maximum cosmic 
ray energy would be E^ax < 2 x 10^^ Zr]iijL33{Wbiowout/Q-7y^^ eV, which is too low, for Z ~ 1. 

Furthermore, ordinary supernova remnants in galaxies appear to have sizes typically 
less than ~ 50 pc (see Matonick and Fesen 1997 and references therein). Therefore, from 
Chevalier's energy-size relation, the energy per event would have to be ~ iQ^os ^^^^ _ jj^dced, 
the standard result, that supernovae create explosions with energy ~ 10^^ ergs, comes from 
precisely these considerations. Requiring the rotational energy per magnetar to be lO^"*^ ergs 
or less requires initial angular velocities no greater than 1.4 x 10^ J^^^^^ sec~^, corresponding 
to a maximum cosmic ray energy not greater than 7 x lO^^Zrji/ 145 eV, which is insufficient 
to explain the highest energy cosmic rays {if Z — 1 — 2). 
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Therefore, the most plausible conclusion is that gravitational radiation must be the 
recipient of most of the initial rotational energy. The ultra-strong magnetic fields in these 
stars make this a viable possibility. From (34), the energy delivered would then he E = 
(l/2)/fi2lnA(fi ^ 0)Wuov.out,GR = 0.75 x lO^^ifiy h^ej) {\n A /2 A) {WMon^out/ 0-7) ergs. Re- 
quiring the interstellar shells created by these electromagnetic explosions to all have radii 
less than 100 pc demands that the neutron stars have equatorial ellipticities e > 2 x 10^^; 
requiring the interstellar remnants of the explosions to have sizes less than 50 pc demands 
e > 4 X 10~^. Prom expression (19), such ellipticities would be expected for interior magnetic 
fields > (6 — 8) X 10^^ Gauss, large but within the realm of recent theoretical suggestions. 
Then more than 90% of the initial rotational energy of the fast rotators radiates away as 
gravitational waves, which still leaves more than enough energy for the UHECR. 

In passing, it is interesting to note that Matonick and Fesen (1997) found a small 
number of unusually large (diameters exceeding 100 pc) supernova remnants in their sample, 
with velocities too large to be readily explained as the consequence of several ordinary but 
coeval supernova explosions. Such a species of supernova remnant clearly is a candidate for 
identification with the electromagnetically driven explosions discussed here. However, given 
the selection effects, it is unclear whether the number of these unusual events is consistent 
with electromagnetically driven explosions from magnetars (or any other exotic explanation) 
as their origin. 



6. Gravitational Radiation 

Since this model requires significant infiuence on the magnetars' spindown by gravita- 
tional radiation losses, with release of most of the initial rotational energy going into this 
so far undetected carrier, the model makes reasonably definite predictions for the detection 
of almost coherent signals by LIGO, VIRGO and other gravity wave detectors optimized 
for khz signals. The dimensionless strain of the narrow band oscillatory signal, measured at 
distance D from the source and observed for a time longer than the coherence time of the 
oscillator Tcoh, so that n 0,Tcoh oscillations can be coherently counted, is (Thorne 1997, 
Brady et al. 1998) 

K = 2.9-^v/^^; (61) 

for simplicity, 1 have averaged over all orientations of the rotation axis with respect to the 
observer, and I have included only the wave modes at the star's rotation frequency, which 
exist for magnetically distorted stars (Bonazzola and Gourgoulhon 1996). The coherence 
time is Tcoh — ^/\^\- 

Since gravitational radiation controls the early spindown [tem ^ tgw)-, Tcoh = 4rGVK[l + 
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{t/TGw)\ and Q = + t/TowY'^ for t < 2Tg ^ l.SIf^el/ /i^^ hours. Then 

^° \ raw ) 



7-3/8 

- 2-5x10--—^^; (62) 

here D20 = -D/20 Mpc. For i > 2rg (rotation periods exceeding 4 msec), the electromagnetic 
torques take over and 

7-2 ^ , r9/4 

/^f ^) = 6.5 X 10-^^^^^y?h^ = 2.9 X 10-^^^^^. (63) 

These strain levels are potentially observable, although the problems of finding such 
quasi-coherent signals in blind searches, with no advance notice of the period, are quite 
formidable. Nevertheless, since the site of the cosmic ray acceleration is at small enough 
radii (r < i?^ ~ 1 AU), the UHE cosmic rays from an individual outburst should arrive 
before the gravity waves fade from view - conceivably, a UHE cosmic ray event, if observed 
with sufficient collecting area to identify a burst of particles, might serve as a marker for 
searches of gravity wave data for quasi- coherent signals with periods in the one to few 
millisecond regime. 



7. Individual Events 



Charged particles of energy in excess of 10^*^ eV must come from nearby For such 
ultra-high energy particles, the GZK losses limit the observable volume to size 
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> 50 Mpc, for maximum energies appropriate to the model discussed here. See expression 
(4) and the immediately preceeding discussion for the values of the characteristic lengths and 
energies. The number of source events per unit time which we can detect, in principle, in the 
ultra-high energy particle spectrum then is Nuhe = {4:71 /3)D^ngi'm"'^^^ ■ This rate is ~ eight 
thousand per year at = 6 x 10^^ eV {De ~ 2.1 Gpc), drops to 800/year at 4 x 10-*^^ eV 
{De ^ I Gpc), but then declines rapidly, at a rate oc E~^ to a few times per decade above 
3 X 10^° eV, assuming the fiducial event rate per unit volume UgU^"'^^ = 2 X 10-^ Mpc-3 
yr~^. The particle fiuence per event is 
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here Dioo = -D/lOO Mpc. / ~ 2 — 3 is a logarithmic function of E and D. 

The model proposed here requires the particles to be accelerated and emitted in an 
equatorial slice with opening angle 2i, with i the angle between the rotation axis and the 
magnetars magnetic moment, into the surrounding Universe - the particles are physically 
beamed. The fraction of the sky into which the particles are beamed is 6 = sin i. If i is 
uniformly distrubuted between and 7r/2, the average beaming fraction is {b) = 2/7i = 0.64. 
The individual fiuence/event would be higher by a factor 1/6. Of course, the rate at which 
events are observed is lower by a factor b, which leads to the average spectrum shown in 
Figure 1 being independent of b. 

Occasionally (once in l/um""^^^ ~ 10^ years), such events must occur in our own galaxy, 
which would create a primary particle fiuence at the earth of order 10^ km~^, for an average 
distance of 5 kpc. Such events might leave tracers in the geochemical record, a topic beyond 
the scope of this study. 



8. Discussion: Related Models 

Blasi et al. (2000) presented the model most closely related to the theory outlined here. 
They were the first to identify the relativistic wind of a magnetized rotator as a potential 
site of UHE cosmic ray acceleration. They assumed the accelerators to be initially rapidly 
rotating pulsars in our own galaxy (fij ~ 3000 s^^), and made the same assumption of the 
charge lost per second, / = qNi = nfl'^/c — c$. Since they assumed normal (if rapidly 
rotating) pulsars, they produced particles with energy ~ 100 EeV only by assuming the 
high energy ions to be fully stripped iron. They made the same construction of the 
particle spectrum as in §3.1 above, for the pure electromagnetic spindown case. They did 
not include possible gravitational wave losses. Since the pulsar spin down time for their 
assumed parameters is long compared to the expansion time of the supernova envelope, the 
pulsar's rotational energy loss does not substantially affect the dynamics of the supernova, 
and most of the UHE particles form at times late enough to avoid substantial radiation and 
stopping losses in the supernova envelope. They gave no consideration to the disruption of 
the envelope due to the pulsars' magnetic pressure on the expelled envelopes. They did not 
use supershells and the sizes of supernova remnants, to constrain the event rate. The large 
charge/particle assumed allowed them to argue that scattering in the galactic magnetic field 
might allow this hypothesized source to be in reasonable accord with the observed UHECR 
isotropy. 

However, UHECR observations support light nuclei, rather than Fe, as the particles 
responsible for the UHE showers {e.g., Bird et al. 1994, Watson 2002). X-ray observations 
of pulsars suggest the atmospheres from which the ions would be extracted are primarily H 
or He, as discussed in §2, not Fe. The energies of protons or a-particles would be too low to 
explain the super-GZK particles in their model. Also, the Larmor radii would be too large. 
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to allow isotropization in the galactic magnetic field. Finally, within the scheme actually 
outlined by Blasi et ai, with each pulsar's magnetic field contained within a closed shell of 
supernova ejecta, the supposed UHE cosmic rays suffer catastropic adiabatic losses. Because 
of the relatively small rotational energy per pulsar assumed in their model, and the slow 
release of that energy {tem ~ 1 year for their parameters), disruption of the shell and escape 
of the wind and its high energy particles takes too long for the UHE particles to survive, 
with ordinary (if rapidly rotating) pulsars as the drivers. 

The metagalactic magnetar model outlined in this paper avoids all of these difficulties, 
as well as adding consideration of the specific mechanism of acceleration in the wind and of 
gravitational wave losses on the rotation history of the pulsar, which can have observable 
effects on the UHE particle spectrum - see Figure 1. 

The theory of supernova envelope disruption and relativistic expansion of a Magnetar 
Wind outlined in §§4.1 - 4.2 have close kinship, and owe much to, the Gamma Ray Burst 
models of Wheeler et al. (2000) and of Inoue et al. (2002). The principal difference from 
the Wheeler et al. scenario is in the disruption of the supernova envelope by Rayleigh- 
Taylor instabilities, rather than assuming a jet punches through the envelope. Inoue et al. 
asummed a rather lower energy rotator, similar to Blasi et al. 's fiducial objects, so that 
Rayleigh- Taylor disruption of the ejecta happens at rather later times. They also assumed 
the Pulsar Wind Nebula in their model would be radiatively efficient so as to construct a 
Gamma Ray Burst model based on nonthermal emission from the Pulsar Wind Nebula, even 
though they also assumed the outflow to be in the form of an electromagnetically driven 
relativistic blast wave, an outlow which may be an inefficient radiator, at least until it 
decelerates to nonrelativistic velocities. In the present model, while there may be an initial 
flash of Compton scattering as the electromagnetic flelds disrupt the supernova shell, the 
subsequent flow is highly magnetized. Such flows have weak shock dissipation (Kennel and 
Coroniti 1984), and may have weak magnetic dissipation (Lyubarsky and Kirk 2001), so the 
conversion of outflow energy to efficiently radiating accelerated particles is likely to be weak. 
However, other processes in a strongly magnetized medium associated with intense electric 
current flow may allow conversion of magnetic energy to emergent photons from the wind 
with unknown efficiency - see Blandford (2002) for comments on some of the possibilities. 

The considerations of laminar surf-riding acceleration outlined in §4.3 have their roots 
in the work of Buckley (1977) and of Contopoulos and Kazanas (2002). The relation of these 
accelerating, force-free wind results to MHD analyses of relativistic winds will be discussed 
elsewhere. 

Magnetars behaving as outlined in this paper probably cannot co-exist with long lived 
{tdisk 3> tgw) fall back disks, as are invoked in collapsar and supranova models of GRBs - 
such structures would probably suppress or greatly modify the flow and acceleration of the ion 
return currents in the rotational equator. However, if that disk is itself strongly magnetized, 
in principle it could play the same role as a UHECR accelerator. The open magnetic flux 
in the magnetar model, ^ = Rl^ ~ 10^^ Maxwells, is similar to that expected in many 
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models of Gamma Ray Burst sources' disks, either due to binary mergers or to fallback 
in a supernova (Narayan et al. 1992, Usov 1994, Paczynski 1998, Vietri and Stella 1999, 
Macfadyen and Woosley 1999). Many of these models involve magnetic fields with magnitude 
> 10^^ Gauss rotating at Kcplerian angular velocity at r > GRschwarschUd = 179M/10 Mq 
km. These objects would have \E' ~ 10^° Mx and $ ~ 10^"^'^ Volts, clearly making them 
available as candidates for UHECR acceleration. Time dependent decay of the disk might 
play a role similar to a magnetar's spindown, in creating a broad energy spectrum of high 
energy particles. Such a model has the advantage that the > petagauss magnetic fields in the 
disks, if both the disks and the strong fields exist, are surely rapidly rotating. The magnetar 
model has the advantage that neutron stars with petagauss surface fields do exist, but we 
do not have independent evidence on how many are born rapidly rotating. 



9. Conclusions 

The model's basic results and predictions are shown in Figure 1. The upper cutoff of 
the spectrum, at £■ = 10'^^ -^Zrji 1133 fll eV, is hkely to be reduced by dissipation of the wind 
as it first breaks free of the magnetar's natal supernova. The estimates of §4.1 suggest that 
escape from the supernova consumes a fraction 1 — Wuowout ~ 0.3 — 0.4 of the magnetar's 
intial energy loss, in turn suggesting that ~ 80% of the highest voltages are available for 
cosmic ray acceleration. If so, then the highest energies in the spectra shown in Figure 1 
should be be taken seriously. 

Agreement with the well determined UHE cosmic ray flux below 5 x 10^^ eV constrains 
the combination Wgeomng2^m4:/ f^ss to be approximately 0.02-0.06, where /iss is the magnetic 
moment in units of 10^^ cgs, z/m4 is the magnetar birth rate in units of 10~^ yr~^, = 
n galaxy/ 0-0'^ Mpc~^ and Wgeom = 0.5 is the inefficiency factor due to the fact that only 50% of 
the magnetic geometries are appropriate to the ion current being emitted into the rotational 
equator, where radiation and adiabatic losses on the UHEGR would be weak. The small 
value of this normahzation required to get a fit, appropriate to the "no GR" and "moderate" 
GR cases considered in Figure 1, implies that the birth rate of magnetars with voltages high 
enough to create UHECR is 5-10% of the overall magnetar birth rate, a conclusion consistent 
with not overpopulating our and other galaxies with too many supershells and overly large 
supernova remnants. 

The high energy ions suffer negligible adiabatic and radiation losses in escaping the 
magnetars' winds and their surrounding, relativistically expanding nebulae. The existing 
observations require voltages $ > 2 x lO'^^/Zrji Volts, or fi^sQ^ > O.OG/rjiZ {Pmitiai < 
2.6(Z?7iyLi33)^/^ ms), with rji = r]/0.1 the fraction of the voltage actually sampled by each 
ion. I suggested in §4.3 that electromagnetic surf-riding in the relativistically strong elec- 
tromagnetic waves in the wind, generated by the oblique rotator, is responsible for the ion 
acceleration, occuring at or within the radius (<S 0.1 AU) where the wave structures no 
longer are frozen into the expanding plasma. That acceleration site is consistent with small 
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radiative losses. 

This model works only if the ions travel from their source galaxies on more or less 
straight lines. That requires intergalactic magnetic turbulence to be small amplitude, or 
the intergalactic magnetic field to be weak, or both - 6B/B < 10~^(10~^ Ga,uss/ ZBigmY^^ 
(§3.4). Such weak scattering introduces negligible time delays into the particles' transport 
from source to observer. 

If the current data are interpreted as providing evidence for the conventional GZK cutoff, 
i.e. the AGASA results are disregarded (Bahcall and Waxman 2003) or are renormalized 
to bring them into accord with the Hi- Res results (Abu-Zayyad et al. 2002b), agreement 
of this model with existing observations requires rather strong losses of rotational energy 
due to gravitational radiation. The "strong OR" curve in Figure (1) shows the effect of the 
rapid spindown depopulating the highest energy end of the injection spectrum. The model 
can replicate the high energy end of the spectrum observed by AGASA only if gravitational 
wave losses are small. The metagalactic magnetar model has no trouble with a spectrum 
that compromises between the highest energy AGASA results and the Hi-Res (without the 
1995 Bird et al. event) results - including that data point favors such a compromise. 

Even when one assumes a birth rate of fast magnetars ~ 5-10% of the overall magnetar 
birth rate, the requirement that the model not overproduce supershells and large supernova 
remnants in our and other galaxies' interstellar media suggests gravitational wave losses are 
an important drain of rapidly rotating magnetars' initial rotational energy. Therefore, one 
expects bursts of almost coherent gravitational waves with millisecond periods and strains 
with magnitude > 10"^-*^ if observed in an experiment designed to find almost coherent 
signals, lasting for several hours and overlapping the arrival of the higher energy particles. 

The model suggests the UHE cosmic rays come from sources whose distribution should 
mimic that of luminous baryons. More specifically, the sources should follow luminous matter 
with stars that are progenitors of core collapse SNc, thus should anticorrelate with large 
galaxy clusters. Events with energy above 10^° eV come from small distances {D < 50 
Mpc), which may allow some imprint of galaxies' large scale structure to appear on the 
isotropy of such events. Starburst galaxies are obvious suspects for especially luminous 
UHECR sources. 

In conclusion, I have shown that bare magnetars in normal galaxies provide a possible 
source for the origin of UHECR, both in total flux and in form of the spectrum. I have 
also given plausible arguments for how the particles accelerated by a magnetar embedded 
in its supernova and its magnetar wind nebula can escape without catastrophic losses, to 
contribute to the particle spectrum observed at the Earth. 

This subject will remain observationally driven, most prominently by the AGASA {e.g. 
Takeda et al 1998), Hi- Res {e.g. Sokolosky 1998, Abu-Zayyad et al 2002b) and Auger {e.g. 
Boratov 1997, Cronin 2001a,b) experiments, which provide a window into the spectrum at 
energies above lO^'^-^ eV. The theory described here suggests that there will be something 
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to see at these super-GZK energies. Also, Auger, with several thousand km^ collecting 
area and with fluorescence and particle detectors recording events simultaneously at the 
same site, will have the opportunity to resolve individual UHE particle bursts, with enough 
statistics to measure a spectrum for one event, at least up to 10^° eV - see §7 for bookeeping 
on this model's predictions. On the theoretical front, the theory of how the acceleration 
actually occurs, in the model presented here and in other schemes, urgently needs substantial 
improvements. 
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A. Appendix - R-mode Gravitational Wave Losses 

Gravitational wave losses modify the accelerated spectrum by shortening the time each 
star spends at a given rotation frequency, thus reducing the flux at energy E{Q). I used grav- 
itational wave spindown rates owing to permanent nonaxisymmetric quadupolar distortions 
of the stars in order to arrive at the spectra shown in Figure 1, and argued that these could 
be generated by strong internal magnetic flelds. The R-mode instability (Andersson 1998) is 
another possible source of gravitational wave emission from a rapidly rotating neutron star. 
With a dimcnsionless mode amplitude a, the spindown rate due to gravitational radiation 
is (Owen et al. 1998) 

-^Rmo^e = [-^ j a, (Al) 

where R is the stellar radius. The same calculation as in §3, with R mode losses replacing 
those due to a static quadrupole moment, leads to 



d'j 




dt \ ^^dn _9 c^I 
4 Ze/i'j 



, 2 

In) """^ dE 



1 + 



(A2) 



replacing expression (23), assuming a = constant during the spindown, as in the model of 
Owen et al.. In (A2), the characteristic energy above which gravitational wave losses affect 
the spectrum now is 
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with Rio — R/10 km. If a were large, Eg{a) might be low enough to yield a steepened high 
energy injection spectrum and a GZK pileup. 

However, recently Arras et al. (2003) provided the theory for a as a function of Q, 
showing that damping in the neutron star's interior, through coupling between R modes 
and internal g modes, keeps a small. In the notation used here, their theory provides 
a = 0.025^4^^. The gravitational wave spindown rate is now proportional to fi^^, in- 
stead of (static quadupole) or (Rossby modes with constant a.) The injection 
spectrum still is given by (A2), but now with (j/^gY replaced by (7/7^)^^^, and with 
Eg = m,ic?^ig = 10^^ Zriifil^/^ / {l45Rio)'^^^ eV. Gravitational wave losses due to the R mode 
instability therefore have negligible effect on the ability of the magnetar model to account for 
existing UHECR observations. They do suggest a very hard cutoff of the spectrum, above 
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